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A 




g 

^gf 

A 

n 


w 

A* 

s 


SgCt) 


Area of a structural member cross section (sq. in.) 

Tributary area of a structure for purpose of computing load 
(sq. ft) 

Area of steel beam cross section, composite beam construction only 
(sq. in.) 

Net area of the rear face of a rectangular structure with openings 
( sq. ft) 

Area of concrete in cross section (sq. in.) 

Area of flange of structural steel beam or girder ( sq. in. ) 

Net area of the exterior front wall of a rectangular structure 
with openings (sq. ft) 

Cross section area of concrete column (sq. in.) 

Gross area of the exterior front face of a rectangular structure 
(sq. ft) 

Area of each portion of the subdivided face of a rectangular struc- 
ture with openings (sq. ft) 

Area of openings in the front wall of a rectangular struct\ire 
(sq. ft) 

Net area of the back wall of a rectangular structure with openings 
(sq. ft) 

Area of reinforcing steel bars parallel to steel beams in composite 
construction ( sq. in. ) 

Area of tension steel in reinforced concrete member (sq. in.) 


Area of one stirrup of the web reinforcement of a reinforced .con- 
crete member ( sq. in. ) 

Area of web of steel beam (sq. in.) 

Area of compression steel in reinforced concrete member (sq. in.) 

Clear distance between flanges of a steel beam (in.)' 

Depth of compression area of concrete (in.) 

Short side dimension of two-way slabs (ft) 

Column spacing in flat slab construction (ft) 

Acceleration of mass (ft/sec^) 

J_T_ 

a a ,, a Acceleration of the g floor mass at times t, 

g,2» g,n-r g,n:fl .j. respectively 

2^ n-1^ rw-l' 

4" V> 

Acceleration of the g floor mass as a function of time 


. a 


n-1" 


a , 
n^ 


nt-l 


Acceleration at time t^, '*'n-l^ ^n' ^n+1' 

respectively 


VII 


EM 1110-345-^13 
1 July 59 


Symbols 


a’ 

a" 

B 



B 

g 

b 


b 

s 

b' 



Distance from edge "a" to centroid of the load on area "A" 

Distance from edge "a" to the centroid of the inertial force 
Peak value of externally applied load 

Elastic stiffness coefficient for soil vtnder rocking foundation 

Equivalent peak value of the externally applied load 

Peak value of the effective load on the g floor mass of a 
multi -story building 

"th 

Peak value of the external load, f (t), on the g floor mass 

s 

Width of flange of steel beam ( in. ) 

Width of conpression flange of reinforced concrete beam ( in. ) 

Long side dimension of two-way slabs (ft) 

Flange width of channel section (in.) 

Width of longitudinal web stiffener for steel beams (in.) 

Width of stem of reinforced concrete tee beam (in.) 

Width of flange of steel beam used in composite construction (in.) 
Distance from edge "b" to centroid of load on area "B" 

Distance from edge "b" to the centroid of the inertial force 

Conpression force developed in reinforced concrete member (kips) 
Static compression column influence factor for shear walls 

Crater depth factor for underground burst 

Dynamic compression column influence factor for shear walls 
Average drag coefficient 

Local drag coefficient 

Maximum ACI Code moment coefficient obtainable for any point on 
a slab (equation 6.52) 

Ratio of maximum resistance to peak load, Cj^ = 

Parameter defined by equation 5*94 
Seismic wave velocity (fps) 

Ratio of load duration to natural period of oscillation, C^ = T/T^ 

Approximate coefficient for determining elastic truss deflection 
(varies with type of truss) ^ 

Ratio of maximum work done to absolute maximum work done, C^zz — 

P 

Local drag coefficient for a cylindrical segment of infinite length 

Value of determinate formed by omitting the row and column con- 
taining the coefficient C. . from the matrix [cj 

1 J 
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D.L.F 


D.I.F 


Distance from neutral axis to the extreme fiber (in.) 

Cohesive strength of the soil as determined by conventional 
laboratory testing procediires (psi) 

The clear distance between compression flange and stiffener 
Velocity of sound in undist-urbed air (fps) 

Velocity of sound in the region of reflected overpressure (fps) 

Diameter of a cylindrical segment or a spherical dome (ft) 

Diameter of the smallest longitudinal reinforcing bar (in.) 
Diameter of spiral core of a concrete column 

Dynamic Load Factor = 

Dynamic Increase Factor for strength of materials 

Depth of burst of bomb below ground stirface (ft) 

Height of biirst of bomb above ground surface (ft) 

Depth from the con^jression face of a beam or slab to the centroid 
of longitudinal tensile reinforcement (in.) 

Least lateral dimension of rectangular column (in.) 

Total depth of steel beam (in.) 

Over-all depth of truss (ft) 

Diameter or width of column capital in flat slab construction (ft; 
Diameter of circular reinforced concrete column (in.) 

Distance between centroids of compression and tensile steel in 
doubly reinforced concrete member 

Modulus of elasticity, including Yoxing’s Modulus for steel and 
soil, of any structural material except concrete 
Energy absorbed by the eq.uivalent system 
Explosive factor for cratering 

Conpressive modulus of elasticity of concrete, kips/sq^. in. 
Initial tangent modulus of elasticity for concrete 
Explosive factor for underground burst 
Maximxjm allowable energy absorption for the g story 
Base of natural logarithms (2.718) 

Eccentricity of load, distance from gravity axis to point of 
application of load (in.) 

Center to center spacing of two trusses 

Strain at failure (^) 

Strain at initiation of strain hardening (^) 

Strain at yield (^) 

Horizontal force 

S\jm of all external loads above the g story, f‘g(‘t) + ^g+i( 4 ) +. 

. fjt) 

Design lateral load on frame at time t^ 
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Symbols 


Fq Summatioii of all external horizontal forces applied to structure, 

including foxindation reactions 

F Pressure factor for underground b\jrst 

Normal component of total passive resistance force (kips) 

F(t) Load applied to a structural element or system as function of time 

f Fiber stress (psi) 

Column critical buckling stress (psi) 

Average stress at design level due to axial load (psi) 

fj Dynamic yield strength for structural steel in composite beam 

construction only (psi) 

f^ Elastic buckling stress due to bending (psi) 

^bd Maximum stress due to bending at design level for steel members 
subjected to eocial load and bending (psi) 

f , Dynamic yield strength of steel (psi) 

'til 

f^(t) Effective load on the g floor of a multi-story structiore as a 
function of time 

f^ Initial static stress (psi) 

f^y Lower static stress (psi) 

f Dynamic yield strength for reinforcing steel for con^osite beam 

construction only (psi) 

f^ Tensile strength of concrete (psi) 

f(t) Load applied to a structural element of system as a function of 

time 

f Upper static yield point stress (psi) 

uy 

fy. Static yield point stress 

f ^ , f„, f_ Factors used in the design of tee beams 
1 ^ 

f^(t), fp(t), f (t) Load applied to the first, second, and g floors 
® of a multi-story building as a function of time 

f^ Static ultimate compressive strength of concrete 

fl Dynamic ultimate compressive strength of concrete 

dc 

G Solidity ratio of truss 

Shear modulus of elasticity (psi) 

G.Z. Ground zero (point on ground below point of detonation of air 
burst) 

g Acceleration of gravity 

H Impulse per unit area. Impulse per foot of length, or total impulse 

(kip-sec) 

Depth of footing or wall below the surface of the r^^ound 
Height of shear wall to center of roof beam 
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Equivalent impulse acting on the equivalent system (kip -sec) 

Total Impulse of the external load on the floor mass (kip-sec) 

Total in^ulse of the external load (kip -sec) 

Height (ft) 

Story height (ft) 

Column height (ft) 

Height of a structirre above groiand (ft) 

Jfeximum thickness of channel flange (in.) 

Clear height of colximn (ft) 

Height of the g"*^^ story of a multi-story building (ft) 

Clearing dimension for each portion of the subdivided front face 
of rectangular structure with openings (ft) 

Clearing hei^t of a closed rectangular structure (ft) 

Vertical distance from roof to point under consideration for buried 

structure (ft) 

Weighted average buildup height of the exterior rear face of a 
rectangular structure with openings (ft) 

Weighted average clearing height for the exterior front face of 
a rectangular structure with openings (ft) 

Weighted average clearing height for the interior rear face of a 
rectangular structure with the openings (ft) 

Weighted average buildup height for the interior front face of a 
rectangular structure with openings (ft) 

Moment of inertia (in.^ or ft^) 


Inertial force of "A" portion of two-way slab 
Average of the gross and transformed moments of inertia 
(in.^ or ft^) 

Inertial force of "B" portion of two-way slab ^ 

Moment of inertia of the gross section ^in. or ft ) 


X raumeiib kj ± xxa^j. o \ ^ 

I^ Moment of inertia of the transformed section (in. or ft ) 

I^ Moment of inertia of rectangular section in concrete tee beam 

design (in.^) / 

Ig Moment of inertia of tee section in concrete tee beam design \^in. j 

I = l[t ] Velocity in5)ulse at time, t^ 

I Ntess moment of inertia of structure about axis of rotation "O' 

j Ratio of distance between centroid of compression and centroid 

of tension to the depth, d 

K Modulus of strain hardening 
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Symbols 


(KE)a 

(KE)e 






K 


Kinetic energy 

Kinetic energy of the actual system 
Kinetic energy of the equivalent system 

Either load, mass, or load-mass factor for slab fixed on fovir sides 

fi = 1, 2, 3 ,...; j = 1, 2, 3>*-*) Stiffness influence coefficient 
(kips/ft) 

Load factor 

Load-mass factor 

ftess factor 

Normal component of passive pressure coefficient accomating for 
the cohesive effect of any soil having an internal friction equal 
to 0 

Normal component of passive pressure coefficient for any soil with 
internal friction angle equal to ^ and with zero wall friction 
developed 

Normal component of passive pressaire coefficient for any soil with 
internal friction angle equal to .0 and with maximum wall friction 
developed 

Resistance factor 

Either load, mass, or load-mass factor for slab simply supported 
on four sides 


KT Kiloton, 1000 tons 

K^ Either load, mass, or load-mass factor for special edge conditions 

Kq Ratio of the maximum average overpressure to the reflected over- 

pressure existing on an inclined roof 

K' Beam equivalent length coefficient used in the design of steel 

beam coluitos 

K" Column length factor used in the design of steel col\amns ariH beam 

columns 

k Spring constant, force required to cause unit deflection of spring 

(kips/ft) 

Soil pressure factor (psi) 


k 

ep 

k-j , k_, 


Soil constant for underground explosion (psi) 

Effective spring constant (kips/ft) 

Equivalent spring constant (kips/ft) 

Spring constant in the elasto-plastic range 

. ..k . .. Spring constants for the first, second, and the g^ 
® stories of a multi-story building (kips/ft) 


Simplified notation of k 
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"th j ... . bh 

-Spring constant of the coupling spring between the g and the 1 
floors of a multi-story building (kips/ft) 

k'. = - k . 
gi gi 

k ,.k k Spring constants of the springs connecting with 

gl -the masses 

Span length of beam or truss (ft) 

Unsupported length of beam (ft) 

Length of- a rectangular structure in the direction of propagation 

of the blast wave (ft) . . 

Length of a cylindrical segment along the cylinder axis ^ft; 

Length of shear wall center to center of column steel 
Spacing of columns in each direction 

Distance from the outside of the front face to the inside rear 
face of a rectangular structure with openings (ft) 

Distance from the front face of a rectangular structure to a point 
under consideration on the roof or sides in the direction of pro- 
pagation of the blast wave (ft) 

Length between sections of zero and maximum moment being considered 
(composite beams) 

Bending moment applied to a section 

Moment of forces on piles about their centroidal axis 
Resisting moment of soil on the footing per unit width 

Bending moment at center line of a beam or slab 

Negative resisting moment in column per foot, flat slabs (kip-ft) 
Positive resisting moment in column strip per foot, flat slabs 
(kip-ft) 

Maximum design moment in a member under axial load, 

Dynamic bending moment 
Elastic dynamic buckling moment 
Maximum bending moment 

Component in a plane perpendicular to edge "a" of the total re- 
sisting moment along the fracture lines bounding area A, two- 
way slabs (kip-ft) 

Component in a plane perpendicular to edge "b," of the total re- 
sisting moment along the fractiire lines bounding area B, two-way 
slabs (kip-ft) 

Negative resisting moment in middle strip per foot, flat slabs 
(kip-ft) 

Positive resisting moment in middle strip per foot, flat slabs 
(kip-ft) 
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M Negative bending moment at support 

Moment of all external forces on structure about axis of rotation 
O „Q„ 

M* Summation of external moments about the point of rotation ex- 

° eluding footing projection moments 



Mp Plastic resisting moment under bending only 

Negative plastic resisting moment in column strip per foot, flat 
slabs (kip-ft) 

M Positive plastic resisting moment in column strip per foot, flat 

slabs (kip-ft) 

Component of the total plastic bending moment capacity along the 
fracture line boundary of area "A" which is in a plane perpendic- 
ular to edge "a," or the total positive plastic bending moment 
capacity for a section parallel to edge "a," two-way slabs (kip-ft; 

M_ Component of the total plastic bending moment capacity along the 

fracture line boundary of area "B" which is in a plane perpendic- 
ular to edge "b," or the total positive plastic bending moment 
capacity for a section parallel to edge "b," two-way slabs (kip/ft; 


“pm 

V> 

Vs 

«Ps 

“Psa 

«Psb 


plastic resisting moment at centerline of beam or slab 
Negative plastic resisting moment in middle strip per foot, flat 
slabs (kip-ft) 

Positive plastic resisting moment in middle strip per foot, flat 
slabs (kip-ft) 

Maximum positive bending moment (kip-ft) 

Plastic resisting moment at support (kip-ft) 

Total negative plastic bending moment capacity along edge "a," 
two-way slabs (kip-ft) 

Total negative plastic bending moment capacity along edge b, 
two-way slabs (kip-ft) 

Yield resisting moment (kip-ft) 

Plastic resisting moment per unit of width of slab (kip-ft /ft) 

Plastic positive bending moment capacity per unit width for short 
span, two-way slabs (kip-ft/ft) 

Plastic negative bending moment capacity per unit width at center 
of edge "a" for long span, two-way slabs (kip-ft/ft) 

Plastic negative bending moment capacity per unit width at center 
of edge "b" for short span, two-way slabs (kip-ft/ft) 
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m 






^2 * 
N 


N 


n 


P 


av 


‘back 



Moment at the intersection of the two-linear portions of the P vs M 
cirrve of a steel heam section (kip-ft) 

Theoretical plastic resisting moment (kip-ft) 

Mass per unit length (kip-sec /ft ) 

Point mass (kip-sec' -/ft) 

In concrete design: 

Total moving mass of stinicture and earth included between footings 
Number of fundamental dimens ioml qiiantities 

Mass of the equivalent system ^kip-sec /ft j 

Mass of structiire considered to rotate as well as translate 

Total mass of the element or structural system under consideration 

(kip-sec^/ft) ^ V 

Total mass of the equivalent system ^kip-sec /ft] 

.m Mass of the first, second, and floors of a multi-story 
® building 


Number of non-dimensional parameters 

Total number of shear connectors required between the points ol 
zero and maximum moment of a composite beam 
Weighted number of piles in group 

Ratio of the length of fixed-edge perimeter to the total perimeter 
Ratio of the length of simply-supported perimeter to the total 
perimeter 


Nimiber of columns in a story j. . • + 

Ratio of modulus of elasticity of steel to modulus of elasticity 

of concrete 

Number of dimensional variables 
Number of stories in a multi-story frame 


Load or force (kips) 

Total dynamic load on slab (kips) 

Panel influence factor for shear walls 

Total load on "A" portion of two-way slab (kips) 

Average axial load acting on each of several columns of a frame 
(kips) 

Local overpressure on the back face of a buried rectangular 
structure (psi) 

Total load on "B" portion of two-way slab (kips) 

Compression mode overpressure (psi) 

Uniform compressive pressure applied radially on arch 
Uniform radial pressure that produces buckling of arch 
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^cyl 


^i) 


'dome 


■ ed 
Pe(t) 

^front 


■ n 


■•P 

^refl 


■ sb 

3 

■ side 

^so 

^soi 

p 

stag 


Local overpressure normal to the exterior surface of a cylindrical 
segment (psi) 

Maximum axial load on column with given M^^ (kips) 

Reduced maximum axial load for long columns 

Deflection mode overpressure (psi) 

Total dynamic axial load 

Local overpressure normal to the exterior face of a spherical 
dome (psi) 

Equivalent concentrated load for equivalent system (kips) 

Elastic dynamic buckling load (kips) 

Equivalent load on an element as a function of time 

Average value of P(t) for the far (or leeward) side of the arch 

llcal overpressure on the front face of a buried rectangular 
struct-ure (psi) 

Peak underground overpressure resulting at a given loeatlon from 
an underground burst (psi) 

Total vertical load (blast plus static) on column 
Force acting at any time, t^ 

Average value of P(t) for the near (or windward) side of the arch 
(psi) 

Dynamic plastic axial load capacity of column 

Reflected shock wave overpressure for angle of incidence of zero 
degrees (psi) 

Reflected shock wave overpressui'e for angle of incidence other 
than zero degrees (psi) 

overpressure existing in the incident shock »ave for any value of 

Ove^rLsure existing in the incident shook wave when t-t^ = t^ 
(psi) 

Local overpressure on the exterior side walls of a rectangular 
structure 


Initial peak incident overpressure (psi) 

ra rt-p c^hock wave formed in the interior of a 

Peak overpressure of the snocK wave 

structure with openings (psi) 
stagnation overpressure; the overpressure 

wSS^the moving air has been brought completely to rest (psi) 
Axial load on the columns of a frame due to the vertical live and 
dead loads (kips) 
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(z? ] The time average of the total vertical load on the column in the 
' time interval "tge ^ ’’^gm 

P Average overpressure on the exterior hack wall of a rectangular 

structure (psi) 

P Average net overpressure acting on back wall of a rectangular 

b-net g-tructure (psi) 

(f \ Peak value of the average overpressure on the exterior back 

( backjmax ^ rectangular structure (psi) 

p Average overpressure on the closed ends of a cylindrical segment 

(psi) 

p Average net overpressure acting on front wall of a rectangular 

f-net g-tructure (psi) 

P Average overpressure on the exterior of the front wall of a 

"front rectangular structure (psi) 

P. Average overpressure on the interior front wall of a rectangular 

structure with openings (psi) 

p Reflected shock wave average overpressure in the interior of a 

i-refl rectangular structure with openings (psi) 

P Average overpressure on the interior roof of a rectangular struc 

i-roof with openings (psi) 

P. Average overpressure on the interior side walls of a rectangular 

i-side structure with openings (psi) 
p Net average horizontal overpressure exerted on a rectangular 

structure (psi) 

p Average net overpressure acting on the roof of a rectangular 

r-net structure (psi) 

p Average overpressure on the front slope of a gable roof (psi) 

'■““f AveSge oveiiressure on the exterior of a reotangnlar structure 

(psi) 

p» Peak average overpressure on the front slope of the first gable 

of a multi-gabled roof (psi) 

P Average net overpressure acting inward on sidewall of rectangular 

s-net structure (psi) 

P . Average overpresaxire on the exterior side walls of a rectangular 

structure (psi) 

P Axial load determined by the intersection of the two linear 

^ portions of the P vs M curve of a steel beam section (kips) 

P(t) Actual load on a structural element as a function of time 
p Ratio of tensile reinforcement in reinforced concrete members to 

concrete area, 

Uniformly-distributed load intensity 
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Symbols 


Q 


*^0 

R 


R, 


R 

R 


dc 


R 


R 


S 


R 


gi 


R 


gm 


R. 

I 

R 


go 

gt 


R 


m 


R 


'mA 


Critical steel ratio for reinforced concrete member 
Equivalent static load for an arch (psi) 

Ratio of volume of spiral reinforcement to the volume of the 
concrete core (out-to-out of spirals) of a spirally reinforced 
concrete column 

Ratio of compressive reinforcement in beams to the concrete area, 
A^/bd 

Statical moment of the section about the centroidal axis 
Strength of a shear connector in composite construction 

Radiant energy on a unit area ^cal/cm j 

Unit drag pressure produced by the incident shock wave (psi) 
Dynamic design bearing pressure on soil ^kips/ft"^] 

Maximum bearing pressure on soil ^kips/ft”^^ 

Maximum unit drag pressure produced by the incident shock wave 
(psi) 

Radius of a spherical dome 

Total resistance of structural element or structural system (kips) 
Dosage of gamma radiation (roentgens) 

Crater radius (ft) 

Total resistance of "A" portion of two-way slabs (kips) 

Crater radius (ft) 

Horizontal static load resistance of shear wall at first cracking 

Horizontal dynamic load resistance of shear wall at first cracking 

Reynolds number of the high velocity wind in the incident shock 
wave 

Resistance of the equivalent system (kips) 

Sing)lified notation of R the resistance function developed 

gjig--*- 

between the g^^ and (g-1)^^ floors of a multi -story building 
The average value of R in the time interval from 0 to t' 

§ O , 

Eesistance function developed between the g and the i floors 
of a multi -story building (kips) 

Maximum resistance of the g'*'^ story of a multi-story building 
(kips) 

Resistance developed in the spring which connects m to the ground 

th ^ 

Total resistance acting on the g floor mass (kips) 

Maximum resistance developed by a structural system (kips) 

Plastic resistance of the "A" portion of two-way slabs where the 
edge is fixed (kips) 
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R. 


mB 


R. 


me 

mf 


n 


R. 


u 


R, 


R 


Im 

ImA. 


R(t) 

Rg(t) 

cv*o 



R 


Plastic resistance of the "B" portion of two-way slabs where the 
edge is fixed (kips) 

Maximum resistance of the equivalent system (kips) 

Fictitious maximum resistance (kips) 

Resistance of a structural element or system at time, t^ (kips) 
Horizontal ultimate load resistance of a shear wall (kips) 

Yield resistance of the structure (kips) 

Resistance in the elastic range (kips) 

Maximum resistance in the elastic range (kips) 

Maximum total resistance in the elastic range - A portion of 
two-way slabs (kips) 

Maximum total resistance in the elastic range - "B" portion of 
two-way slabs (kips) 

Resistance as a function of time 

Resistance of the g story columns as a function of time 

Time variation of the resistance of the equivalent single- 
degree -of -freedom dynamic system for the g™ story 

Dosage of nuclear radiation without shielding (roentgens) 

Resistance as a function of displacement 

Radius of gyration of section (in.) 

Resistance per unit length of a beam or per unit area of slab 
Ratio of web reinforcement » Ay^bs 

Distance from point of explosion (ft) 

Roentgens 

Ratio of steel reinforcing placed perpendicular to the steel beam 
( in composite construction) in excess of that required to carry 
the slab bending stresses 

Average value of the resistance of the 8 ., story when the relative 
displacement between the g"^^ and the g-l''"“ story is negative 


S Section modulus 

SE Strain energy absorption (kip-ft) 

(SE)^ Strain energy absorption of actual structure (kip-ft) 

(SE)^ Stia,in energy absorption of equivalent structure (kip-ft) 

S' Section modulus about the weak axis 

s Coordinate spcis at 45 ° angle to x axis in x-y plane 

Spacing of stirrups and spacing of ties in reinforced concrete 
members (in.) 
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Symbols 


Distance along s axis to centroid of loading (ft) 

Distance along s axis to centroid of inertia force (ft) 

Duration of the external load (sec) 

Resultant tension force (kips) 

Fundamental period of vibration for complete spherical shell of 
centroid thickness (sec) 

Natural period of oscillation or fundamental period (sec) 

Lowest natural period of circular arch with pinned or fixed ends 

'til 

Natural period of the i nrade of oscillation (sec) 

Time or rise of a step-load ( sec) 

T Time duration of the external loads on the first, second. 


IX^ 2X 


and g floors of a multi-story building (sec) 

T py . . »T y Time duration when the effective loads on the first, 
second, and g"*-^ stories are negative (sec) 


Period of oscillation of a fictitious system representing the 
g'th story (sec) 

Thickness of concrete slabs (in.) 

Thickness of flange of beams (in.) 

Time measured after the arrival of the incident shock wave (sec) 
Time variable 

Thickness of deep beam web 

Time of arrival, or time required for the shock wave to travel 
from the point of explosion to the chosen location (sec) 

Average flange thickness of standard steel beam (in.) 

Time required for overpressure on the rear face of a closed 
rectangular structxrre to rise from zero to its maximum value (sec) 

Time required to clear the front face of a structure from the 
reflection effects (sec) 

Time displacement factor; the time required for the shock front 
to travel from the frontmost element of a structure to the point 
or surface under consideration (sec) 

Time at which the limiting elastic deflection is reached 
Flange thickness of a WF steel beam (in.) 

Time at which the limiting elastic deflection of the g*^ story 
is reached 

Time required for the vortex generated at the front face of a 
structure to travel a distance L' across the structure (sec) 

Time required for maximum displacement of element or structure 
to occiir (sec) 

Duration of the positive phase of the incident shock wave (sec) 
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Time of rise; time required for the overpressure in the incident 
shock wave to rise from zero to its maximum value (sec) 

Thickness of stiffener (in.) 

Time required for the shock to pass over the length of the 
structure (sec) 

Time of rise of an undergro-und oveirressure pulse (sec) 

Web thickness of a steel beam or channel (in.) 

Time required to reach yield point of material (sec) 




Time sequence 


*lm^ ■^2m' 


n* n+l 

. .t Time at which the maximum absolute displacements of the 
first, second and g'*'^ floor masses are reached 




At = t 
n 


Time required for the overpressure on the exterior rear face of 
a rectangular structure with openings to rise from zero to its 
maximum value (sec) 

Time required to clear the front face of a structure with openings 
from reflection effects (sec) 

Time at which the relative displacements and respectively, 
are equal to zero 

Time immediately before t^ 

Time immediately after t^ 

Time lag 

Time interval used in numerical analysis 

- t A time interval 
n+l n 

Velocity of the incident shock front (fps) 

Velocity of the shock front of the shock wave formed in the 
interior of a rectangular structure with openings (fps) 

Particle velocity in the incident shock wave (fps) 

Bond stress per unit of surface area of bar in reinforced con- 
crete design (psi) 

Dynamic reaction (kips) 

Total shear (kips) 

Total vertical force on piles 

Total dynamic reaction along one edge "a," two-way slabs (kips) 
Average vertical shear in length, L’ (kips) 

Total dynamic reaction along one edge "b," two-way slabs (kips) 
Total column load in flat slab design (kips) 

Maximum shear capacity of web of deep beam section 
Maximum vertical shear force (kips) 
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Symbols 


Vortex velocity (fps) 

Shear stress (psi) 

Dynamic shear yield strength 

7 ft \ Velocity of the g floor at time, t 


= V ^t \ Velocity of the g fl 
.n g^ 

Maximum shearing stress (psi) 
= Velocity at time, t^ 


v+ Velocity at tf 

n n 

V£ Velocity at t^ 

Horizontal velocity of axis of rotation " 0 " 

v^t^^ Velocity of mass at time t^ 

V Static shear yield strength of steel (psi) 

^ Average unit shear stress in concrete (psi) 

v(t) Velocity as a function of time 

W Total energy yield of an atomic bomb expressed KT of TNT 

required for an equivalent total energy yield 
Weight (lbs) 

Total load on element of structure (kips) 

Work done (ft-lbs) 

Dynamic peak blast load 

W Work done on the actual system (ft-kips) 

Si 

W Work done on the equivalent system (ft-kips) 

^ "til 

W Maximum work done on the g floor mass (ft-kips) 

O » 

W Absolute maximum work done on the g floor mass (ft-lbs) 

W 

m Maximum work done on the equivalent system by the equivalent load 

Wp Fictitious maximum work done on the equivalent system 

W(t) Work done as a function of time 

w Uniformly distributed load (kips/ft) 

Length of channel shear connector ( in. ) 

Width of front face of rectangular building (ft) 

w^ Weighting factor for vertical and batter piles 

X Relative displacement in a story of a multi-story building (ft) 

^gm Maximum relative displacement in the g^^ story (ft) 

^gs Relative displacement of the mass g when the average external 

loads are applied statically 

^1? Xg, X Relative displacement in the first, second, and g^^ stories 
of a multi-story building (ft), X = x - x 
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"til 

X Limiting elastic relative displacement for the g story column 

(ft) 

X Distance the incident shock wave travels after impinging on the 

frontmost element of a cylindrical segment or spherical dome to 
any point being considered 

Deflection of a structural element or system (ft) 

Distance of any pile from the oentroidal axis 

X Collapse deflection 

Llmitir^ elastic deflection (ft) 

X Forced solution of a dynamic system for a given external load 

f til 

X Absolute displacement of the g floor at time t 

go 

X Maxinium absolute displacement of the g floor (ft) 

gP^ th 

X . X . X , Absolute displacements of the g floor at time, 
g,n-l* g,n^ g,tt+l ^ t , and t , 

n-1' n^ iH-1 

X Maximum displacement (ft) 

m 

x^ x^, x^^ Displacement at time t^_^, t^, and t^^, respectively 

X Initial displacement (ft) 

° Horizontal displacement of axis of rotation "0" 

X Displacement of an elastic system subjected to the peak load B 

® acting statically (ft) 


x^tgj Displacement of a mass at time t^ 


£ , X , X Absolute displacement of the first, second, and g floors 
^ ^ ® of a multi -story building (ft) 

ic Distance to centroid parallel to x axis (ft) 

ic Horizontal velocity of axis of rotation "0" 

o 

id Horizontal acceleration of axis of rotation "0" 

id Acceleration in the x direction (ft/ sec^ j 

y Displacement or deflection of a structural system (ft) 

y Deflection of actual element (ft) 

£L 

y Midspan deflection of actual element (ft) 

y^ Midspan deflection (ft) 

y Limiting elastic deflection (ft) 

® Deflection of equivalent system (ft) 

y Limiting deflection in the elasto-plastic range (ft) 

•'ep 

y^ Maximum displacement (ft) 

y Displacement at time, t^ (ft) 
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Symbols 




a 


n 


u 


Vertical distance from axis of rotation ”0'* to mass centroid of 

rotating mass, m . of structure 
r 

Deflection at midspan of column strip - flat slab (ft) 

Deflection at midspan of middle strip - flat slab (ft) 

Maximum elastic range deflection (ft) 

Maximum elasto-plastic range deflection (ft) 

Distance to the centroid parallel to y axis (ft) 

Vertical distance from axis of rotation "0" to mass centroid of 
total moving mass, m 

Acceleration in the y direction (ft/sec^) 

Velocity in the y direction (ft/sec) 

Velocity of displacement of the center of the beam (ft/sec) 
Increment of displacement (ft) 

Plastic modulus of the cross section (in.^) 

Plastic modulus of the cross section in the weak direction (in.^ 


Depth of the resultant compressive force in concrete tee beam 


Angle of incidence between the normal to the surface and the 
direction of propagation of the blast wave (degrees) 

Central angle of the arch 

Deflection coefficient for two-way slabs 

Design load ductility reduction factor 


Angular acceleration of structure about axis of rotation 
(radlans/=e=2) 


Dimensionless ratio = 0.5P /l^-7 
Ductility ratio 


" 0 " 


Maximum bending deflection due to the application of a \anit load 
(ft/kip) 


Shear wall lateral deflection at first cracking 


Sub-8trea clearing factor for rectangular struct\ares with openings 

Maximum shear deflection due to the application of a unit load 
(ft/kip) 

Shear wall lateral deflection at ultimate resistance 
Deflection at theoretical yield (ft) 

Angle of inclination to the horizontal of a gabled roof (degrees) 
Parameter used to define the location of a point on the surface 
of a cylindrical segment or a spherical dome (degrees) 

F.nd rotation of structural element 
Angular displacement 
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Joint rotation at bottom of column 
Rotation of beam at midspan (radians) 

Angular displacement of structvire about axis of rotation "0" 
Joint rotation at top of coliamn 

Viscosity of air in the incident shock wave ^ Ib-sec/ft ) 
Coefficient of friction 

Poisson's ratio , 2 4\ 

Mass density of air in the incident shock wave ( lb-sec /in. j 

/ 2/ 4 2/ 4\ 

Mass density of soil ^lb-sec /in. , or lb-sec /ft J 

Critical buckling stress (psi) 

Parameter used to define the location of a point on the surface 
of a spherical dome 

Phase angle in the transient response (radians) 

Internal friction angle (degrees) 

Phase angle of the i^^ mode (radians) 

Unit weight of soil (ibs/cu. ft) 

Angular velocity of structure about axis of rotation "0" 
(radians/ sec ) 


XXV 



EM 1110-3i»-5-^13 
1 Jtily 59 


T^nuals - Corps of Engineers 
U. S. Army 


ENGINEERING AND DESIGN 

DESIGN OF STRUCTURES TO RESIST THE EFFECTS OF ATOMIC WEAPONS 

WEAPONS EFFECTS DATA 

lETOODUCTION 


3 _Oi FURPQPIR AHD SCOPE. This marcoal is one in a series Issued for the 
gpTt Hp nr»(=> of engineers engaged in the design of permanent type military 
structures required to resist the effects of atomic weapons. It is appli- 
cable to all Corps of Engineers activities and installations responsible 
for the design of military construction. 

fphft material is based on the resiilts of full scale atomic tests and 
analytical studies. The problem of designing structures to resist the ef- 
fects of atomic weapons is new and the methods of solution are still in the 
development stage. Continuing studies are in progress and supplemental 
material will be published as it is developed. 

The methods anrl procedixres were developed throu^ the collaboration 
of many consultants and specialists. Much of the basic einalytical work was 
done by the engineering firm of Ammann and Whitney, Ifew York City, xmder 
contract with the Chief of Engineers. The Massachusetts Institute of Tech- 
j^ology was responsible,^ under another contract with the Chief of Engineers, 
for the conipilation of material and for the further study and development 
of design methods and procedures. 

It is requested that any errors and deficiencies noted and any sug- 
gestions for Improvement be transmitted to the Office of the Chief of Engi- 
neers, Department of the Army, Attention; ENGEB. 

3-02 EEEEEEIICES. Manuals - Corps of Engineers - Engineering and Design, 
containing interrelated subject matter are listed as follows: 

DESIGN OF STRUCTURES TO RESIST TEE EFFECTS 
OF ATOMIC WEAPOHB 

EM 1110-3lt-5-4l3 Weapons Effects Data 

EM mn -RkR-klk Strength of Materials and Structural Elements 
EM 1110- 345 -4i 5 Principles of Dynamic Analysis and Design 
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EM 1110-3ii-5-^l6 
EM lllO-3i)-5-^l7 
EM 1110-3^5-^18 
EM 1110-3^5-^19 
EM 1110-3i)-5-^20 
EM 1110-3^5-^21 


Structural Elements Subjected to fiy irn t ml c Loads 

Single -Story Frame Buildings 

Multi-Story Frame Buildings 

Shear Wall Structiires 

Arches and Domes 

Buried and Semiburied Structures 


I 


a. References to Material in Other Manuals of This Series. In the 
text of this Ttiamiai references are made to paragraphs, figures, equations, 
and tables in the other manuals of this series in accordance vith the 


number designations as they appear 

in these manuals. The first part of the 

designation which precedes either s 

i dash, or a 

decimal point. 

identifies a 

particular manual 

in the series as 

shown in the table following. 


EM 

paragraph 

figure 

equation 

table. 

1110-3ii-5-^13 

3- 

3. 

< 

;3. ) 

1 

3 . 

1110-3l<-5-^l^ 

4- 

4. 

{ 


1 

4. 

1110-3i<-5-^15 

5- 

5. 

( 

!5. ] 

1 

5 . 

1110-3l<.5-Ul6 

6- 

6. 

< 

' 6 . ] 

» 

6. 

1110-31+5-^17 

7- 

7. 

1 

I- ! 

I 

7 . 

IIIO-3I+5-I+I8 

8- 

8. 

1 


1 

8. 

IIIO-3I+5-I+19 

9- 

9. 

1 

(9. I 

) 

9. 

1110-345-420 

10- 

10. 

10. 

I 

10. 

1110-31+5-1+21 

11- 

11. 

(11. : 

) 

11. 


b. Bibliography. A bibliography is given at the end of the text. 
Items in the bibliography are referenced in the text by numbers inclosed Ixi 


brackets. 

c. List of Symbols. Definitions of the symbols used throu^out this 
man'ual series are given in a list following the table of contents. 

3-03 RESC!ISSI0ir> (Draft) EM 1110-3^5-^13 (Part XXIII - The Design of 
Structures to Resist the Effects of Atomic Weapons, Chapter 3 ” Weapons 
Effects Data. ) 

3 - 0 I GENERAL. Atomic weapons may be exploded in the air, at ground level, 
or at various depths below the ground svirface. The types of protected 
structures whose design are covered in this series of manuals are surface, 
shallow covered, and shallow buried structures. The major effects con- 
sidered in protection against atomic weapons are; 

( 1 ) Air blast 

( 2 ) Nuclear radiation 

(3) Thermal radiation 

(4) Cratering 

( 5 ) Ground shock. 
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^ Underground "burst and cratering cause insignificant effects on a 

structure xmless the structure is in close proximity to the hurst. These 
effects have been given very limited coverage in this manual. 

The air burst is the major factor in the determination of the re- 
quired strength of a blast-resistant structure. To evaluate the loads on 
a structure it is necessary to develop the critical blast loading pattern, 
relative to each exposure of the surface of the structtire. Pressure in- 
tensities will vary with respect to time, reckoned from the Instant of 
arrival of the shock front. 

Nuclear radiation is an important consideration in the design of 
protective structures. In some designs radiation shielding may govern the 
required thickness of concrete and earth cover whereas the blast effects 
will be the determining factor relative to the designs for structural 
strength and stability. Consideration should be given to both direct and 
residual radiation. Incendiary properties and heat insulation from thermal 
radiation should be given attention in selecting suitable materials for 
construction of protective structures. 

^ 3-05 BLAST WAVE PHEINOMEHA.. a. Infinite Atmosphere. To simplify the 

presentation of air blast phenomena, the explosion in an atmosphere of in- 
finite extent is first considered. Almost immediately after the detonation 
occxirs, the expansion of the hot gases initiates a pressure wave in the 
surrounding air as represented roughly by the curve in figure 3*l(8')* This 
shows the general nature of the variation of the overpressure (pressure 
above atmospheric) with distance from the esplosion at a given instant 
during the early stages of shock formation. As the pressxire wave moves 
outward from the center of the explosion, the following, or inner, part 
moves through a region which has been previously compressed and heated by 
the leading, or outer, parts of the wave. The pressure wave moves with the 
velocity of soxuid, and since this velocity increases with the temperature 
and pressure of the air throu^ which the wave is moving, the inner part of 
the wave moves more rapidly and gradually overtakes the outer part, as 
shown in figure 3.1(b). The pressure wave front thus gets steeper and 
steeper, and within a very short period the pressure change becomes abrupt 
forming a shock front, as indicated in figure 3*l(c). At the shock front 
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there Is a sudden increase of pressure from | 
normal atmospheric to the peak shock pressure 
and the advancing shock front he comes a moving 
wall of highly compressed air. 

Initially, in the hot central region of 
the explosion the pressure exceeds atmospheric 
hy a factor of many hundred thousand. The 
pressure distribution behind the shock front 
is somewhat as illustrated in figure 3«l(c)* 

It shows the peak overpressure at the shock 
front, dropping rapidly in a relatively small 
distsuace to a value about one -half the shock 
front overpressure. The magnitude of the 
overpressure is essentially uniform in the 
central region of the explosion. 

As the expansion proceeds, the pressure 
distribution in the region behind the shock 
front gradually changes. The overpressure is no longer constant but drops 
off continao\isly nearer the center. At later times when the shock front 
has progressed some distance from the center, a drop in pressure takes 
place below the initial atmospheric value at the center and a suction phase 
develops. The front of the blast wave weakens as it progresses outward 
and its velocity drops toward the velocity of sound in the undisturbed 
atmosphere. The sequence of events Just described is depicted in figure 
3.2. This shows the overpressure dis- 
tribution in the blast wave as a func- 
tion of the distance from the e3q)losion 
at different stages in the expansion. 

When the suction phase is well de- 
veloped, the overpressure in the blast 
wave resembles the heavily drawn curve 
in figure 3.2. The pressure variations 
in the blast wave from this time on can 
be considered in two different ways as 

4 



Figure 3.2. Variation of overpressure with, 
distance at various times 



(c) 


Figure 3.1. Overpressure during 
early stages of shock 
front formation 


Blost Wavt<< 


NtgofiyePhose 


Positivt Phase 
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shown in figures 3*3 and 3 . 1 +a. 

First, the heavy curve of figure . 

3 1— Uo 

3.2 has "been redrawn in figure 3-3 to S shock Fronf^^,^^ 

show the variation of overpressvLre | “ 

with distance at a given time. The r^';tW7;;nr^|r.itiv. pho,, | 

symbol represents the peak 

overpresstire, or shock intensity, in 

pounds per square inch and repre- Figure 3.3. Variation of overpressure with 

sents the velocity of the shock front ® 

in feet per second. The arrows vmder the curve show the direction of move- 
ment of the air mass, or blast wind, in the positive and negative phases. 

^ .j The peak overpressure in the positive 

I p V phase is higher than the maximum 

I \. overpressure in the negative phase. 

' Consequently, the blast wind is of 

Timooftorootonotion higher velocity in the positive phase 

Vnrtnfinn nf^,,^rr,r^ssurewith than in the negative, or suction. 


Oittonct from Centtr of Explosion 

Figure 3.3. Variation of overpressure with 
distance at a given time 


Time offer Detonation 


Figure 3.4a, Variation of overpressure with 
time at a given location 


phase. 


Second, the pressuire variations in the same wave may be considered 
at a point located a given distance from the essplosion as indicated in 
figure 3 .i+a, which shows the variation of overpressure with time at a fixed 
location. The symbol t is the time of arrival, or the time in seconds 
for the shock front to travel from the center of the explosion to the given 
location, t is the duration in seconds of the positive phase, and p 
is as previously defined. 

The rate of decay of the positive phase varies with respect to the 

magnitude of the peak overpiassure. In figure 3*^^ "the rate of decays 

for various magnitudes of peak overpressures is expressed in terms of 

the ratio P/P to the ratio t/t , where P^ is the overpressure at 
'so o' s 

time t after the arrival of the shock front. The peak negative over- 
pressure is approximately one-eighth of the peak positive overpressure 
and the duration of the negative phase is approximately four times the 
duration of the positive phase of* the blast. 

The peak overpressure P in free air before the shock front 
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Figure 3.5. Peak overpressure for 1-KT free air burst by 
2W assumption (taken from reference[3^) 
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Radius of 
Shock Front 


strikes the ground Is a function of the distance from the point of hurst 
and the yield of the weapon. Values for 1 -KT bomb given in figure 3-5 are 
from, curves given in reference [36]. which are based on the supposition that 
a contact surface burst of W kilotons energy is equivalent in blast 
characteristics to an e3q)losion of 2¥ kilotons hi^ in the air and prior 
to reflection. 

b. Finite Height of Burst. The discussion in the preceding para- 
graph deals with the air blast from an atomic bomb e3!p)loded in an atmos- 
phere of infinite extent. A discussion follows relative to the influence 
of the surface of the earth on propagation and attenuation of the air blast. 

If the bomb is detonated at a hei^t h above the surface of the 
earth, the shock front will have the general configuration indicated in 

figure 3*6 during the brief time inter- 
^^N^center of val before it impinges upon the surface., 

/ Explosion 

Radius of / \ The shock front is the same as if propa- 

Shock Front jr\ vl^ 1 

I f gated in an atmosphere of infinite ex- 

\ / ^ 

Shock tent* 

Front N. X 0 

^ ^ •• A short time later the radius of 

Earth's Surfoce---^^ ® 

the shock front becomes greater than 

_ - , r.T . If 1 . 13^ Sind a portion of the shock front 

Figure 3.6. Shock front before striking 

ground surface impinges upon the earth ' s surface and 

is reflected back, forming the reflected shock front illustrated in 
figure 3 * 7 * The arrows at the incident and reflected shock fronts indicate 
the direction in which the blast waves 
are traveling. The location of the re- 

fleeted shock front is roughly deter- Sho'e*" / \ 

mined by drawing an arc with center lo- \ 

cated a distance h below the earth's I Reflected 

surface and directly under the point of A yL 

detonation, Joining the points of inter- surface 

section of the incident shock front with ^ » nr ^ r 

tigure 3. 7 . bhock front reflection at 

the earth’s surface. The symbol a ground surface, a less than 45^ ± 

represents the angle of incidence of the shock front with the earth’s sur- 


Eorth's Surface- 


Figure 3.6. Shock front before striking 
ground surface 


Incident 
Shock 
Front / 


Reflected 
Shock 
Front y 


\ / 


/Earth's 

77P7//7?7/7^>/7//7Mw//m 


-y Surface 


Figure 3. 7. Shock front reflection at 
ground surface, a less than 45 ° ± 
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is a function of the 
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incident shock overpressure P aud the angle of Incidence 

SO 


a [20]. 

The reflected shock front in figure 3.7 travels through the atmos- 
phere at a hi^er velocity than the incident shock because it is traveling 
through a region of greater than atmospheric pressure and gradually over- 
takes merges with it to form a single shock front called the Mach Stem 

as shown in figure 3.8. The fused shock front thus fomed is normal to and 
tmvels parallel to the earth's 
surface. The Mach Stem forma- 
tion is initiated when the 
angle of incidence o: of 
the shock front becomes greater 
than approximately 45 degrees. 

Once formed, the height of the 
Mach Stem gradually increases 
as the radius of the shock 
front becomes greater. 

The region on the earth's 
surface within which a is 

less than 45 degrees, and no Figure 3.8. Shock front reflection when a is 

Mach Stem is present, is called greater than 45°± 

the region of regular reflection, while the region for which a is greater 
than 45 degrees and a Mach Stem is present is called the region of Mach re- 
flection. The importance of the Mach Stem phenomenon is that it results 
from the fusing of two shock fronts to form a single shock front of higher 
overpressure and of greater destructive potential to structures located in 

its path. 

The shock front velocity is a function of the peak overpressure 

P ; its value for standard atmospheric conditions is plotted in figure 3*9- 



so 


t of the shock front at a given location is a func- 
a 


The time of arrival 
tion of the distance from the center of the explosion and the total energy 
yield of the weapon. The duration on the grotind surface of the positive 
phase t of the blast wave is a function of the peak overpressure P^^ 
and the total energy yield of the weapon. Values of the duration of the 
positive phase are plotted in figure 3 • 10a and in figure 3 • dOb . 
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Reflected ovei^ressure ratio plotted in figure 3*11 as 

a function of the angle of incidence o: of the shock front. This figure 
applies to an inclined shock front striking the earth’s surface or to a 
vertical shock front striking a reflecting surface such as a wall of 
structure . 

The peak overpressure in the blast wave (shock front overpres 

sure) at the ground surface plotted against distance from ground zero for 
various homh sizes is given in figure 3.12a for a typical air hurst and in 
figure 3.12h for a surface hurst. The typical air hurst is at a height 
above the ground which will maximize the target area covered by the 12-psi 
blast pressure level. This height probably cannot be attained by manned 

aircraft for the larger size bombs. 

Por other heists of burst may be obtained from charts in 

reference [1+9] , where the peak oveipressure is plotted as a function of the 
height of burst and the horizontal distance from ground zero. 

The variation of the overpressure in the blast wave on the ground 
with time can be detemined in the same manner as for the bomb burst in an 
infinite homogeneous atmosphere. Equation 3.1 plotted as figure 3.1<-‘i is 
valid for ideal wave shapes which occur at approximately 10 psi or lower. 

For higher oveirressuie levels the variation of overpressure with time at 
the ground level Is usually Irregular, due to thermal effects and to other 
modifying Influences of the ground surface. Reference [49] contains data 
on these irregular wave forms. However, due to lack of data and the com- 
plexity of calculations involving irregular wave forms, the ideal wave foims 
that are given in figure 3.^c should be used for practical design purposes. 

The dynamic pressure q due to the motion of the air particles in 
the blast wave is given by the relation q = pu"/2 where p is the mass 
per unit volume of the air and u is the velocity of the air particles. 

The peak values of q at the ground surface plotted against distance from 
ground zero for various bomb sizes are given in figure 3.23b i'be peak 
values of q plotted against the peak overpressure P^^ are given in 

figure 3 . 23 a 
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The rate of decay of the dynamic pressure is a function of the peak 

overpressure. Figure 3.4d shows variation in dynamic pressures for ranges 

of overpressures as a function of 0/0.2^ S'! a fraction of time t/t^j 

where P is the overpressure and q. is the peak dynamic pressure and 
so o 

q is the dynamic pressure at time t , after the arrival of the shock 
front . 

3_06 SCALUTG blast PHENOMEM. It has heen found that air blast phenomena 
such as the pressure and duration at different distances are related, for 
different strength homhs, according to the ratio of the cube root of the 
yield expressed in TNT. These relations are referred to as scaling laws. 

These scaling laws state that if a given peak overpressure is experi- 
enced at distance r^ from an explosion of a bomb of total energy yield 
the same peak oveipressure will be experienced at distance r^ from 
the explosion of a bomb of total energy yield where: 



The same scaling laws also state that while the peak pressures from 
the two bombs are equal to the two radii r^ and r^ , the durations of 
the blast pressure waves at the two points are different. If the duration 
of the positive phase of the pressure wave from the first bomb is t^^ at 
radius r^ the duration of the positive phase t^^ pressure wave 

from the second bomb at distance r^ will be: 


t 


o2 



(3.3) 


LOADING ON RECTANGULAE STBUCTURES WITHOUT OPENINGS 


3_07 loading on structures. The manner in which the blast wave loads a 
structure varies with the ratio of the distance of the structure from 
ground zero to the height of the burst. The loading on a structure located 
within the region of regular reflection (figure 3*13) is greatly different 
in character from that for a structure located beyond this region. Since 
little or no data are available within the region of regular reflection and 
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Reflected 




Figure 3. 13. Structure located in region Figure 3. 14. Structure located in region 

of regular reflection of Mach reflection 

since a surface structure located in this region, would he subject to such a 

hi^ overpressure as to preclude a practical design^ the procedures for the 

determination of loads are restricted to surface structures located outside 

the region of regular reflection (figure 3*1^) where the Mach Stem is high. 

enou^ to cover the structure. 

The problem is that of computing the loading on a stiucture due to 
the Impingement upon it of a blast wave traveling parallel to the surface 
of the earth. The stiucture is considered as being oriented with one face ■ 
normal to the direction of propagation of the shock wave, since such an I 

orientation produces the most severe loading on the structural elements. 

For the design of a structure so located, it is necessary to predict the 
overpressures which would exist on various portions of the structure as a 
function of time t measured after the shock front strikes the front wall. 

The procedures developed in the following paragraphs for the deter- 
mination of the loads on different types of structures are based on experi- 
mental results obtained from TNT and atomic explosions and from shock tube 
data. The computational methods presented are based on an interpretation 
of these experimental results and the equations given in connection with 
the loadings on a structure result from curve- fitting of measured data and. 
should not be regarded as representing the development of any new theory. 
Nevertheless, every attempt has been made to make the empirical relation- 
ships consistent with such theory as now exists. 

The following paragraphs are devoted to a general presentation of tlie 
procedures for the deteimination of the loads on structures of various 
types. Pertinent quantities are introduced and explained, and Justification 
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for the methods employed is given as well as the bases for their concep- 
tion. The general presentation for each type of structure is followed by a 
paragraph giving a step-by-step procedure by which the loads can be com- 
puted. Illustrative numerical examples are presented for the closed 
rectangular structures. 

In using procedures presented herein, inconsistencies in loads com- 
puted may be encountered, particularly in the case of buried and semiburled 
structures and with peak overpressures greater than 25 psi. The lack of 
experimental data and an adequate theory in these regions requires use of 
engineering judgment and some degree of conservatism. Extrapolation and 
use of load data for peak oveipressures greater than, 25 psi are strictly 
an expediency and must be recognized as such by the engineer. 

a. Diffraction and Drag Loading. The loading on aboveground struc- 
tures resulting from the air blast produced by an air or surface burst may 
be considered to consist of a diffraction phase and a drag phase. 

The diffraction phase of the loading is the term given to the initial 
phase of the blast loading on a structure when the reflected pressures as- 
sociated with the air blast are acting on the structure. The time required 
for the blast wave to surround the structure completely and the presence of 
large reflected pressures on the front wall cause net lateral loads to be 
exerted on the structure as a whole in the direction of travel of the blast 
wave. The local and differential forces which act on the structure during 
the initial stages are defined as the diffraction phase loading. 

The drag phase of the loading is the term given to the second phase 
of the loading on a structure due to the mass and velocity of the air 
particles in the blast wave after the envelopment of the structure by the 
wave front and the reflection effects have decayed. This phase of the 
loading is most important on open structural frameworks and on structures 
having small dimensions such as stacks, poles, chimneys, etc. 

In general, the diffraction phase of the loading can be neglected and 
only the drag phase considered if the minimum dimension of the structure 
perpendicular to the direction of travel of the shock wave and the dimen- 
sion of the structure parallel to the direction of travel of the blast 
wave are less than approximately 5 Because of such short dimensions. 
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the duration of the diffraction phase is very short The larger the loaded i 
area the greater is the time required for the transient reflection effects 
to decay and reach the relatively steady condition of the drag phase® The 
longer the duration of the reflection effects the more Important is the 
diffraction phase of the loading. In all cases except when a precursor 
exists^ the pressures on any portion of the stiucture existing during the 
diffraction phase are greater than during the diag phase® 

Large quantities of heated dust may he raised hy the thermal radia- 
tion before the arrival of the blast wave® The air near the ground is also 
heated by the thermal radiation and this may cause a precursor wave to form, 
ahead of the main shock. Test results indicate that when precursors form 
there is usually a dust-laden atmosphere present and the dynamic pressures 
are much larger than when there is no precursor. These high dynamic pres- 
sures are usually observed in the high oveipressure regions. 

3-08 LOADING ON CLOSED RECTMGUIAR STRUCTURES® The type of structure for 
which the loading predictions of this section are applicable is illustrated, 
in figure 3 .15® The behavior of the blast wave upon striking a closed rec- | 
tangular structure is depicted in figure 3*l6(a)^ (b), (c), and (d). This 
figure illustrates the position of the shock front and the behavior of the 


reflected and diffiacted waves at successive times during the passage of 
the blast wave over the center, portion of the structure of figure 3*15® A.s 
the shock front strikes the wall of the building (figure 3 • 16 (a)) a re- 


flected blast wave is formed and 
the overpressure on this wall is 
raised to a value in excess of the 
peak overpressure in the incident 
blast wave. This increased 



Figure 3^ IS, Closed rectangular 
structure 



Figure 3, 16, Behavior of blast wave along ce/z- 
ter portion of closed rectangular structure 
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overpressure is called the reflected overpressure and is a function of the 
peak overpressure in the incident blast wave and the angle of incidence of 
the shock front with the front wall which is zero degree in this case. At 
the instant the reflected shock front is formed, the lower oveipressure 
existing in the incident blast wave adjacent to the top edge of the front 
wall initiates a rarefaction wave (figure 3.16(b)), or a wave of lower 
overpressure than that which exists in the reflected blast wave. This 
rarefaction wave travels with the speed of sound in the reflected blast 
wave toward the bottom of the front wall. Within a short time, called the 
clearing time, the rarefaction wave enfeebles the reflected blast wave to 
such an extent that its effect is no longer felt and reduces the overpres- 
sure existing on the front wall to a value that is in equilibrlxim with the 
hi^-velocity air stream associated with the incident blast wave. The 
overpressure on the front wall when equilibrium with the high-velocity air 
stream is reached is equal to the stagnation oveipressure at the base of 
the front wall and an overpressure somewhat less than that in the blast 
wave at the top edge of the front wall. The stagnation oveipressure is de- 
fined as that overpressure existing in a region in which the moving air has 
been brought completely to rest causing the pressure intensity to be in- 
creased by the loss in momentum of the air particles. 

At some time after the shock front strikes the front wall of the 
structure, equal to the length of the structure divided by the shock front 
velocity, the shock front reaches the back edge of the structure and starts 
spilling down toward the bottom of the back wall (figure 3 •18 (c)). The 
back wall begins to experience Increased pressures as soon as the shock 
front has passed beyond it. The effect is first observed at the top por- * 
tions of the back wall and proceeds toward the bottom. A vortex, which is 
a region of air spinning about an axis at a high speed with low overpres- 
sures existing at its center because of the venturi effect, is created on 
the back wall and grows in size, traveling toward the base from the top 
edge and also moving away from the wall (figure 3 *18(1) and figure 3«17)* 
The maximum back wall oveipressure develops slowly as a result of the vor- 
tex phenomena and the time required for the back wall to be enveloped by 
the blast wave. 
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Figure 3.17. Photograph of a blast wave passing over a square block. 

Light and dark bands represent contours of constant pressure 

As the shock front passes beyond the front wall (figure 3.16(b)) the 
overpressure exerted on the roof of the structure is initially raised to a 
value nearly equal to the overpressure existing in- the incident blast wav6 
However, the air flow caused by the pressure difference between the re- 
flected oveipressure on the front wall and the blast wave overpressure on 
the roof causes the formation of a vortex along the top edge of the front 
wall. The vortex travels with a gradually decreasing intensity along the 
roof of the structure (figure 3*l6(c)) at a slower rate than the shock 
front velocity. It causes a decay of the overpressures built up by the in- 
cident blast wave. After the passage of this vortex, the higher overpres- 
sures in the blast wave again become dominant and cause a second build-up 
of overpressures along the roof. 

If a horizontal section through the structure is considered as shown 
in figure 3 *18, it will be found that the blast wave produces effects simi- 
lar to those produced on the roof, front, and back walls as indicated in 
figure 3 * 16 . The major effects which should be noted are: 
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(1) The formation of rarefac- 
tion waves at each end of the front 
wall of the structure which travel 
toward the center to relieve the 
hi^er reflected pressures on the 
front wall; 

(2) The diffracted wave 
fronts at the rear edges of the side 
walls which travel toward the center 
of the back wall; 

(3) The vortices which fom on 
the side walls near the front edges 
and travel slowly toward the back; 

(4) The vortices which form at 
the ends of the back wall and travel 
slowly toward its center. 

In determining the loads on a 
structure, it is convenient to use 
the instant at which the shock front 


I Incident 
U^Shock Front 


Vortex 


•Shock Front 


End Wall 


End Wall 


Rarefaction 
Wove 


Reflected 
Shock Front 



Shock Front 


iffracted 

Shock 

Front 



(d) 


Figure 3.18. Behavior of blast tvave along horizon- 
tal section through closed rectangular structure 


impinges on the front face as a reference for time (t = 0). In adopting 
this convention, it is necessary to introduce a time -displacement factor 
t which is the time required for the shock front to travel from the front 
face of the structure to the surface or point under consideration. Figure 
3.19 illustrates this convention. In this figure, t - t^ is the time 
after the shock front has passed a given location on the structure and de- 
termines the overpressure in the 

incident blast wave at that location. 
For example, the time-displacement fac- 
for a given location on the roof a 
distance L' from the front face of 
the structure is t^ = • 

The total loading on any face of 
a cubical structure is equal to the 

Figure 3.19. Time-displacement 

factor convention algebraic sum of the respective 


Time Scale 
Shock Front 
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oveipressures or reflected oveipressures combined with the dynamic pres- ' | 
sures. The later pressures are due to the wind effects of suspended parti- 
cles and are referred to as "drag" loading which may result in either posi- 
tive or negative corrections depending on the orientation of the surface 
with respect to the direction of travel of the shock front. 

a. Average Front Wall Oveipressure At the moment the inci- 

dent shock front strlkes'the front wall, the overpressure on the frorrt wall 
is immediately raised from zero to the reflected oveipressure Re- 

flected overpressure is plotted in figure 3 ‘20 for zero angle of incidence 
as a function of the peak overpressure of the incident shock front. The 
value of P from figure 3.20 is the same as P for a = 0° using 

figure 3.11. The incident shock front continues its motion over the top of 
the structure, while the reflected shock front moves away from the front of 
the building in the opposite direction. Initially, the air pressure be- 
tween the reflected shock front and the front wall is higher than the pres- 
sure behind the incident shock front. This causes air to move around to | 
the sides over the top of the structure into the lower pressure zone be- 
hind the incident shock. The rarefaction waves thus formed move from the 
end and top edges toward the center of the front face with the speed of 
sound for the pressure existing in this region of reflected overpressure flj, 
For the range of shock strengths used in the Princeton shock tube ex- 
periments [4] which correspond to peak shock oveipressures of 2 to 50 psi 
in a standard atmosphere, it has been found that the time required to clear 
the front face of reflection effects is determined by the dimensions of bhe 
front face and the peak oveipressure of the incident shock wave. This 
clearing time t is given by the relation: 

*0 " f3.4) 

where 

h' = clearing height, taken as the full height of the front face 
or half its length whichever is the smaller 

Crgfi = velocity of sound in the reflected region, plotted as a func- 
tion of the peak overpressure of the Incident shock wave in 
figure 3.21 

During the time required to clear the front wall of reflection ef- 
fects, the average overpressiire on the front wall decreases from the 
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0 20 40 60 80 100 
Peak Incident Overpressure, llo(psi) 

Figure 3,20, Reflected overpressure vs peak incident overpressure for normal reflection 
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reflected oveipressure to a value given ly the following equation: 




(3.5) 


Where P = overpressure in the incident "blast wave at the front wall at 
® any time t - t^ as given hy the equation 

Equation ( 3 . 6 ) is plotted in figure 3*22, and tabulated in table 3»1> where 
the ratio Pg/Pg^ is given as a function of (t - t^)Ao* This equation is 
similar to equation ( 3 *l) with (t - t^) substituted for t. 

In paragraph 3-05b the dynamic pressures associated with high over- 
pressures were discussed. These pressures were evaluated in figure 3 •23a. 
The rates of overpressure decay and dynamic pressure decay were set forth 
in figures 3.4b and 3.4c. 

Design requirements for aboveground frame structures are usually 
for peak overpressures of 10 psi or less. The use of equations (3*6), 
( 3 . 7 a), and 3 . 7 c following and tables 3.1 and 3.2 will be applicable 

and convenient. 

•The quantity q in equation (3-5) is the dynamic pressure due to the 
motion of air particles in the incident blast wave at any time (t - t^) and 
is given by q = pu 2/2 where p is the mass per unit volume of the air and u 
is the velocity of the air particles. Values of dynamic pressures vs over- 
pressures for the theoretical (without precursor) and precursor condition 
are given in figure 3.23a. It will be noted that the values are about the 
same for each case up to about 10-psi overpressure. The dynamic pressure 
vs distance for miles for various weapon sizes is given in figure 3 . 23 b. 


q = 14.7 


1 (P /i4.7)‘ 


15 


(3.Ta) 


1^1 + Y 

Peak dynamic pressure q for use when overpressures are less than 
10 psi or where there are no precursor effects considered is given by the 
equation P ^ 




1 t i (^30/44.7) 


(3.TI) 
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Table 3.1. Blast Wave Overpressure Ratio vs Time Ratio 
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Figure 3.23b. Dynamic pressure for surface burst 
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To ottain q. as a function of time for overpressures less than 
10 psi it is convenient to use the ratio plotted in figure 3*24 emd 

tabulated in table 3*2 as a function of (t - from the following 

ejcpression 




(3.7c) 



Time Ratio, ( t - t^,) /t^ 


Figure 3»24, Dynamic pressure ratio vs time ratio 
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Table 3.2. Dynamic Pressure Ratio vs Time Ratio 
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0.17 
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0.20 

0.21 
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0.24 

0.25 

0.26 
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0.28 
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0.873 

0.835 

0.798 

0.762 
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0.552 

0.527 

0.503 

0.480 

0.458 

0.437 
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0.379 

0.361 

0.344 
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0.39 
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0.56 
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For the front wall overpressures, t, = 0, hence t - t, = t, where 

'a ' d ^ 

t measures the time after the Instant at which the shock front impinges on 
the front wall. 

Using the above quantities, the curve of average front wall overpres- 
s-ure versus time may be determined. Figure 3.25 shows a typical 

front wall average overpressure curve. The curve a-b-c is defined by 

equation (3- 5)* Point d, the re- 
flected overpressure, is connected by 
a straight line to point b, the aver- 
age overpressure at time t - t^ * t^. 
Since the resulting discontinuity at 
point b is not compatible with 
actml behavior of the loads on the 
front face, these two curves are 
smoothed by fairing in curve e-f as 
shown. The curve of average over- 
pressure versus time on the front wall is then defined as ctirve o-d-e-f-c. 

f^oregoing method for determining the curve of average overpres- 
sure versus time on the front wall is based on the assuugption that the time 
required for the overpressure in the Incident shock wave to rise from zero 
is negligible. 

If this is not the case, the time required for the peak overpressure 
to rise from zero, called the time 
of rise t^, must be taken into 
account. This is accomplished by 
first contputing the average front 
face overpressvire for a zero time 
of rise as described above. Curve 
o-d-e-f-c of figure 3.26 illvis.- 
trates the average front wall 
overpressure computed for a zero 
time of rise. Locate point g 
with time coordinate t^ over- 
pressure coordinate Join points o and g with straight 



Time, t 


Figure 3*26, Graphical procedure for evaluating 
reduced average front wall overpressure due 
to time of rise of incident shock wave 



Figure 3*25, Average front wall overpressure 
vs time— closed rectangular structure 
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intersecting d-e-f at h. The average front wall overpressure for time of 

rise t is given by curve o-h-e-f-c. 

me overpressure on the front wall of a structure is not uniformly 
distributed. The maximum value occurs at the mid-point of the base and the 
value occuTS along the edges. Those portions of the front face 
nearest to the edges are cleared of the reflection effects in a shorter 
time than the remainder of the front face and the overpressure existing at 
those points is lower, following the clearing stage. The net effect of 
this vertical and horizontal variation is negligible and of questionable 
value for design purposes. Hence, the front wall loading is assumed to be 
distributed unifomnly over the front wall surface. For the same reasons, 
the rear wall loading given next is assumed to be uniformly distributed 


over the rear wall. 

b. Average Back Wall Overpressure \ack‘ 
time-displacement factor t^ is • L/U^, where^'L is the length of the 
building in the direction of propagation of the shock and is the 

shock front velocity. 

When the shock front crosses the rear edge of the structure, the foot 
of the shock spills down the back wall. The overpressures on the back wall 
behind this diffracted wave are considerably less than those in the inci- 
dent blast wave due to the vortex which develops at the top and travels 
down the wall. A period of time longer than that required for the travel 
of this diffracted shock to the bottom of the wall must pass before the 
back wall average oveip)resS'ure reaches its peak value. The time required 
for this build-up to occur, t^, measured from the instant at which the 
shockwave reaches the back wall is equal to kh'/c^, where h' is the 
clearing height of the back wall, taken equal to either the full height 
of the back wall or half the width of the building, whichever is the 
smaller, and c^ is the velocity of sound in undisturbed air, 1115 fps. 

The peak value of the average overpressure on the back wall after 
this build-up has been completed is: 

(P J .P^(l/2)ri+(l-3)a'M (3-8) 

'back' max sb ' L -* 

where P is the incident blast wave overpressure at the back wall at 
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on the hack wall which occ\irs at time t»t, + t, ; 6 = 0. 5P /l4.7: and 

d b so' 

e = 2.7183 = base of natural logarithms. It is assTimed that P in the 

s o 

incident blast wave does not diminish in strength as the wave passes over 
the structure. 


Figure 3 •27a illustrates the variation of the ratio P, , /P with 

back' s 

time. For times in excess of t = t^ + t^, the relation is given by the 
following equation: 

^’back _ ^^back^max ^ _ ^^back^max ^ ~ ^^d 

^sb _ ^sb ^o ” \ 

where t^ is the duration of the positive phase. 

Figure 3 •27b shows a typical back wall average overpressure curve. 




For ^ ^bock Lineorly 

from Zero ot t=tj,to the Value Indicoted 
at t*tb +t^. 


Figure 3.27a. Average back wall 
overpressure ratio vs time— 
closed rectangular structure 



Time , t 

Figure 3.27b. Average back wall overpressure 
vs time— closed rectangular structure 



Time, t 

Figure 3.28. Net horizontal overpressure vs 
time— closed rectangular structure 


c. Average Net Horizontal Over- 
press-ure Considering as positive 

all overpressures exerted on the struc- 
ture and directed toward the interior^ 

^net ~ ^front ” ^back (3-10) 


vhere 


^net’ ^front' 


back 


given in terms of time t. Pigtire 3.28 
illustrates graphically the variation 
of the net average horizontal over- 
pressure with time. 
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Area Affected 
by Presence 
of Vortices 


Blest Wave 


d. Local Roof Overpressure procedure developed for pre- 

diction of oveii)ressures on the roof of a structtire is tased primarily on 
curve-fitting of the Princeton and Michigan shock tube data [4 throu^ 13 ]. 
Pressures calculated by this procedure have been found to agree satisfacto- 
rily with overpressure records obtained in the GEIEEIIHOUSE structure testing 
program and the Sandia HE tests [6]. 

During the passage of the blast wave across the structure, low pres- 
sure areas develop on the roof and side walls due to vortex formation as 
indicated in figure 3.I6 and figure 3 - 18 . Shock tube data indicate that 
this vortex detaches itself from the front edges and moves across the 
structure with gradually Increasing 

b g 

speed. Vortices are formed all along the 
front edges of the roof and sides of the 

structure, that is along a-b-c-d of ^ 

figure 3.29. At corners b and c 

where they are aligned at 90 to each d e 

other, the vortices tend to Interfere Figure 3.29. Areas or. roof and sides 

with each other and to move away from most affected by vortex action 

the roof and wall surfaces. This results in a diminution of the effects 
of the vortices at these edges which proceeds along the axes of the vor- 
tices from corners b and c as the vortices move toward the rear of the 
roof. Consequently the regions on the roof and side walls which are 
strongly affected by the vortices do not extend to all edges of the roof, 
but ai^ triangular or trapezoidal in shape as indicated in figure 3 * 29 . 

C.L. -?< ^ These observed results indi- 

, cate a variation of overpressures on 

Bla»tWave ^^zont 2 \ 

\ I the roof in a direction parallel to 

‘ l/4 Zone I Shock front in addition to a 

Ii Zone 2 variation along the roof in the dl- 

J 

rection of propagation of the wave. 

Zont 3 

Althou^ this lateral yariation is 

Length, L— 

a smooth one^ the roof and side 
Fie^re 3.30. Location of loading zones on . ^ 

roof and sides of structure walls have been divided into three 

zones as illustrated in figure 3-30 “to facilitate computation of local 


Figure 3.29. Areas on roof and sides 
most affected by vortex action 


Zone 3^ 

Blast Wove ^^:^ne 2 
Zone I 


: — r~ 

1 L^4 

Zone 1 

Zone Z 



Zone 3 


Length, L 
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pressures. Although only the location of the zones on the roof is dis^ 
cussed helov, the zones on the side vail are similarly- located as shown in 

figure 3*30* 

Zone 1. This zone is a strip on the roof extending from the sides 
toward the ^nter line a distance equal to one-quarter of the length L of 


the building. 



Figure 3.31. Local roof or side wall over- 
pressure ratio vs time, Zone I— closed 

rectangular structure 


for any point in Zone 1 is equal to 
the overpressure-time curve for the 
incident blast wave displaced in 
time by factor t^ « L‘/U^o This 
is consistent with observed vortex z 

CL 

disintegration as illustrated in ^ 

figure 3-29- ? 

o 

Zone 2. This zone is a strip ^ 
on the roof of width equal to one- t 
quarter of the length of the struc- | 

ture measured from the edge of Zone ^ 
1 toward the center of the roof* o 

Figure 3-32 is a plot of the 
ratio a 

distance L* from the front edge 


Figure 3*31 is a plot of 

9 "the ratio of the local 
roor s^ 

overpressure at any point on the 
roof in terms of time t to the 
overpressure in the incident blast 
wave where time- 

displacement factor^ equal to the 
time required for the shock front to 
travel a distance ^ the distance 
from the front edge of the roof to 
the point under consideration. 


Figure 3*31 indicates that the 
local overpressure versus time curve 



Figure 3.32. Local roof or side wall over- 
pressure ratio vs time^ Zone 2— closed 
rectangular structure 
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and 


r „ P incorporate the time- 

roof s 

t is given by equation (3«ll)- 
m 



C.L. 


Note! t(i= lVUo 


of the roof. Here again, P 
displacement t^ = _ 

Zone 3. This zone includes all points on the roof not included in 
Zones 1 and 2. Figure 3-33 is a plot of the ratio 
points in this zone. Zone 3 is 
essentially a region of two- 
dimensional shock phenomena for 
which considerable data are avail- 
able [ 4, 5, 9 , 10, 11, 12, and 13]. 

The effects of the vortex 
travel across the roof of the 
structure are confined to Zones 
2 and 3. The difference in the 
overpressiire existing in these 
zones is due to the assiimed 
severity of the vortex effects. 

As the shock front passes 
over the point being considered, 
the local overpressure is raised 
to the overpressure in the inci- 
dent blast wave. This equality between local and incident blast wave over- 
pressures is maintained until the vortex developed at the front edge of the 
roof detaches itself and moves toward the rear edge, causing a decrease 
from the overpressiire in the incident blast wave at the point being con- 
sidered. The local overpressure reaches its minimum value at the time that 
the vortex is over the point in question. The time required for the vortex 
to travel the distance L' is 

t = L'A (3.11) 

where v is the vortex travel velocity, obtained from a reduction of the 
data in reference c^] and given by the relation 



Figure 3.33. Local roof or side wall over- 
pressure ratio vs time. Zone 3— closed 
rectangular structure 




042 + 0.108 


The value of the ratio of ^roof/^s 


f)u 

L / o 
at time 


t * t 


m 


is 




(3.12) 


(3.13) 
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For points located in Zone 2 the minimum value of eqmtion (3*13) is 
0.5; for points located in Zone 3 the minimum value is zero. 

As soon as the vortex has passed the point in question, the local 
overpressiire starts to build up until at time t » t^^ + 15 h'/UQ (where h' 
is the clearing hei^t of the structure) it is once again eqvial to the 
overpressure in the incident shock wave. 

e. Average Roof Overpressure Having established the varia- 

tion of local overpressure, the average overpressure on the roof at any 
time could be obtained by the summation of the local overpressiire curves 
over the entire roof at a given time t. If there were no lateral varia- 
tion of local overpressures, such a method could be used to develop a 
general procedure applicable to all structures. However, the presence of 
this lateral variation conrpll cates any procedure to such a degree that only 
the limiting case for Zone 3 loading is considered. 'The ratio of the aver- 


P 2 

NOTE* P”* 0.9+ 0.1 (l.0“~); •xc«pt P" moy not oxcood 1.0 
P‘ • 2n-(^ + 1 Xf ) '^or 0.8+ 0.125 (2- ,^) * 


whichtvor is tmolltr; ixcopt P* moy not bo loti 

« than zoro. 

0 . 



Timo, t 

Figure 3.34. Average roof or side wall 
overpressure ratio vs time. Zone 3 — 
closed rectangular structure 


age roof overpressure 

the overpressure of the Incident ! 

blast wave existing at the ■ 

center of the roof P is 1! 

s p 

plotted in figure 3*34 for times 

in excess of time t » l/u 

' o 

where P incorporates the 
s 

time-displacement t^ = 1/211^. 

The average oveip)ressure on the 

roof varies linearly from zero 

at time t » 0 to the value of 

Proof gdven in figure 3*34 at 

time t » assun^tion 

that expressed as 

a proportion of P is valid 
s 

only if the time required for 


the shock front to travel the length of the building is small coiiq)ared to 
the duration of the positive phase. For values of (L/u^)/t^ less than 0.1 
this assunption is reasonable. 

A further restriction which must be iiii)Osed is due to the importance 


42 


3 -08f 


EM 1110 - 3 ^ 5-^13 
1 July 59 


of the lateral variation of overpressures on the roof. If the width of the 
structure normal to the direction of travel of the shock wave is greater 
than twice its length, the average roof overpressures as determined by 
figure 3 . 3 ^ are satisfactory. For a structure whose width is less than its 
length, the average roof overpressure for times greater than time t » l/u^ 
is more correctly given by the relation 

P „ » P ( 3 * 1 ^) 

roof s 

where P^ incorporates the time -displacement t^ * The average 

roof overpressirre varies linearly from time t » 0 to the value calculated 
by equation (3.14) at time t = L/Uq- For structures which are approxi- 
mately square in roof plan, neither method is the more correct and an aver- 
age of the two should be used. If a more accurate value of the average 
overpressure on a portion of the roof, such as a roof slab, is desired, 
this may be obtained by computing the local overpressure at two or three 
points on the slab and obtaining the wei^ted average as a function of 
their common time t. 

f. Local and Average Side Wall Overpressure and On 

the basis of test data it has been observed that the sides of a structure 
are loaded in the same manner as is the roof, and that figures 3.31 to 3.34 
determined for the roof of a structure apply equally well to the sides. 

The side walls are divided into three zones,- as shown in figure 3*30^ 

for the determination of local overpressures Figures 3-31 to 3*33 

present these local overpressures in the form Fgj_( 3 ^g/Fg; where Fg^^g and 

P incorporate the time-displacement t » L'/U • L' is the distance from 
S u. w 

the front edge of the side wall to the point under consideration. 

The average overpressure on the side is given as the ratio ^slde/^s 
in figure 3-34, where P^ is the overpressure in the incident blast wave 
with a time -displacement t^ * L/2U^. This figure is applicable only if 
the height h of the side wall is greater than half the length L. If 
this restriction is not satisfied the average overpressure on the side is 
more correctly given by the relation 

P = P ( 3 . 15 ) 

side s 

3-09 PROCEDURE FOR COMPUTATION OF LOADS PIT CLOSED RECTAtTGULAR STRUCTURES. 
The following paragraphs present a series of step-by-step general procedures 
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for the camputation of the blast- wave overpressirres on the various por- 
tions of a closed rectangular structvire. 

Step 1. Given data: Size of weaponj distance of proposed structure 
from ground zeroj exterior dimensions of structures. 

Step 2. Determine for the given data from figure 3 •12a or 

3.12b using the scaling relations of paragraph 3-06 as necessary. 

Step 3. Determine t^, knowing by applying the scaling rela- 

tion of equation (3.3) to figure 3- 10a or 3.10b. 

Note that all the quantities listed in Step 1 are not necessary to 
determine the loading on any structure if P^^ and t^ plus the exterior 
dimensions of the structiire are given. 

a. Average Front Wall Overpressure vs Time. Step 1. Determine P^ 
for a sequence of times t^, t^^ t^^ etc., from table 3*1 or figure 3*22 


Step 2. Determine knowing P^^, from figure 3.20. 

Step 3. Determine the clearing height of the front wall h’ . This- 
is eq-ual to the full height of the front wall or half its width, whichever 


is the smaller. 


Step k. Determine 
Step 5 • Determine 


Crefl^ Isnowing 


P , from figure 3*21. 
so 


from the relation 


'refl* 


Pg + 0.85q 


Step 6. Determine q^, knowing P^^, from figure 3.23- 
Step 7- Determine q for a sequence of times t^, tg, t^^ etc., 
from table 3.2 or figure 3.2^1-. 

Step 8. Determine (P + 0.85q) for the sequence of times used in 

s 

Steps 1 and 7 above. The following headings systematize the calculations 
of Steps 1 to 8. 

time, t t/t^ q/<j^ g P^ + 0.85q 

(1) (2) (3) W (5) (6) (7) 

Column (l) is a set of times ranging from t » 0 to t » t^ suffi- 
cient in number to allow drawing of the curve of P + 0.85q. 

s 

Column (3) is values obtained from table 3*2 or figure 3*24 for the 
values of t/t^ ih column (2). 

Column (5) is values obtained from table 3.I or figure 3.22 for the 
values of t/t^ in column (2). 
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Step 9' Plot the curve of time as illustrated in figure 

3 . 25 . 

h. Average Back Wall Overpressvire vs Time. Step 1. Determine 
Imcwing PgQ, from figure 3 . 9 . 

Step 2. Determine t^, from the relation t^ » 

Step 3' Determine t^, from the relation t^ » kh^/c^, vhere c^ » 
1115 fps and h' is either the full height of the hack wall or half its 
length, whichever is the smaller. 

Step *i-. Determine P^^ at time \ table 3 .I or 

figure 3 . 22 . 

Step 5 . Determine (P, , ) at time t = t , + t. from the rela- 

oacK max d d 

■ I [1 + (1 - ’*<=« 6 - 

The average overpressure on the rear face varies linearly from zero 


at time t 


■to (Pv 1 ) 

' back'max 


at time t » t^ + t^^- 


Step 6 . Determine the values of Pijack^^s ^ sequence of times 


in excess of t = t^ + t^ from the relation 


back 

P 


(P ) (P ) 

^ back max ^ ^ back max 

P + i - p 

sb sb 


t - (t^ + \) 

"o-S 


Step 7- Determine the values of P^ for the sequence of times in 
Step 6 , and hence the values of Th® following tabttlar headings 

systematize the work in Steps 6 and 7- 

tima, t t - (t - 


back 


( 1 ) ( 2 ) ( 3 ) W ( 5 ) (6) (7) 

Column ( 1 ) is a sequence of times ranging from t = t^ + t^ to 

t = t + t, . 
o d 

Column (4) is a set of values obtained from table 3*1 or figure 3*22 
for the values listed in column ( 3 ). 

Colxmin ( 6 ) is the values calculated from the relation in Step 6 for 
the times given in column (l). 

Step 8 . Plot the curve of vs time as illustrated in figure 

3.27b. 

c. Average Uet Horizontal Overpressirre vs Time. Step 1. Knowing 
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„ ^ and P, , for a sequence of times t^, t_, etc., determine 

front _ back ^ I-" 2 ' 3 ' 

. = - P.,. I for each value of t. 

net front back _ 

Step 2. Plot the cvirve of P^^^ vs time as illustrated in figure 

3 . 28 . 

d. Local Roof Overpressure vs Time. The step-by-step procedure 
listed below is for a point located within Zone 3 (see figure 3 -30) of the 
roof. The procedure for Zone 2 is similar except that figure 3*32 applies 
instead of figure 3 - 33 , and the value of the ratio at time 

t « tm (step k) may not be less than 0 . 5 . The procedure for Zone 1 in- 
volves a simple con 5 )utation of at various times using table 3*1 or 

figure, 3. 22 as illustrated in figure 3 * 31 * 

Step 1. Determine t^ = L' /U^. 

Step 2. Determine the vortex velocity v from the relation 
V » (o.Oi+2 + 0.108 ^ j 
Step 3* Determine t = L'/'v* 

HI 

step 4. Determine the value of the ratio P ^ at time t » t 

from the relation Pj.QQf/Pg = except that the 

minimum permissible value of this ratio is zero. The ratio ^j-oof/^s 

creases linearly from 1.0 at time t = ('t^ + ^roof^^s 

t and increases linearly to 1.0 at time t » t + 15h'/U j for all other 
m HI o 

times the value of the ratio is 1 . 0 . 


Step 5« Determine the value of the ratio ^poof'^^s ^ sequence 

of times between time t » t , and time t » t , + t . 

d do 

Step 6 . Determine P for the sequence of times in Step 5- 

s 

step 7- Determine Pj.qq^ from Steps 5 and 6. The following tabular 
headings systematize the work in evaluating P^qq^* 


time, t t - t . (t - t,)t P /P P P ^P P 

' d ' d' o S' so s roof' s roof 

CD ( 2 ) ( 3 ) CM (5) (6) (T) 

Column ( 1 ) is a sequence of times ranging from t * t^ to 

t a t, + t . 
a 0 

Column (4) contains the values obtained from table 3 .I or figure 
3.22 for the ■values in column ( 3 )* 

Column ( 6 ) con-tains the values obtained in Step 4 above. 
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Step 8. Plot the cxirve of vs time from the values in column 

( 7 ) as illustrated in figure 3 . 33 . 

e. Average Roof Overpressure vs Time. The step-hy-step procedure 
listed below is for a structure for which figure 3 * 3 ^ is applicable^ as 
explained in paragraph 3 - 08 e. 

Step 1. Determine t^ = L/PTJ^. 

Step 2. Determine the value of P „/p at time t = t, + L/2U 
^ roof' s d ' o 


from the relation Pj.QQf./Pg * 0*9 + 0.1 (l.O - 


so 


<1.0 


Step 3- Determine the value of ^pQQf/^g ‘time t » 51>/U^ from 


the equations P^Q^f/^s = 7 \irj 


iL 


or P „/P » 0.5 + 

roof' s 


0.125 (2 - whichever gives the smallest value, but not less than 


zero. 


Step il-. At time t » 5L/U^ + 15h'/U^ ^1^® of ^roof/^s 


one. ^roof'^^s linearly between the values conputed in Steps 2, 

3 , wtiH 4 , and remains equal to one for times greater than t » 51 -/^^ + 

15h'/Uo. 

Step 5 . Determine the values of Pj.QQf/Pg fo^ ® sequence of times 
greater than t = L/U by interpolating between the values conputed above. 
Step 6 . Determine P for the sequence of times used in Step 5* 
step 7 . Determine Pj^qqj from Steps 5 and 6 . The following tabiilar 
headings systematize the work in evaluating P^qq^.- 


t - t . 


(t - t^)/t^ 


P /P 
s' so 


^roof/^s 


roof 

(7) 


time, t 

( 1 ) ( 2 ) ( 3 ) W (5) ( 6 ) 

Column ( 1 ) is a sequence of times ranging from time t = L/U^ to 

time t = t + t , . 

0 d 

Column (4) contains the values obtained from table 3*1 or figure 3-22 
for the values in column ( 3 )* 

Column ( 6 ) contains the values obtained from Step 5 for the times in 
column ( 1 ) . 
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Step 8. Plot the curve of vs time as Illustrated in figure 

3.3J+. Note that varies linearly from zero at time t » 0 to the 

computed value at time t « 

f. Local Side Wall Overpressure vs Time. The local side wall over- 
pressure- for a point in Zone 3 of the side is coii5)uted by use of Steps ! to 
8 of paragraph 3-09d. The procedure for Zone 2 is similar except that 
figure 3.32 applies instead of figure 3-33, and the value of the ratio 

P /P at time t » t (Step 4) must not be less than 0.5. The proce- 
roof' s in. 

dure for Zone 1 involves a simple con5)utation of P^ at various times 
using table 3.I or figure 3-22 as illustrated in figure 3 •31* 

g. Average Side Wall Overpressxare vs Time. The average side wall 
oveip)ressure for a structure for which figure 3-34 is applicable is com- 
puted by use of Steps 1 to 8 of paragraph 3-09© . 

3-10 NUMERICAL EXAMPLE OF COMPUTATIONS OF LOADS ON CLOSED BECTANGUIAR 
STRUCTURE, 



Note! Length of Bldg. Parallel to Shock 

Front i. 200 ft. From figure 3.10a t^^ = 0.262 

Figure 3.35. Exterior proportions— closed (l KT for 10 psl) . 

For an IB-KT weapon, t . 18^/3 
(0.262) = 0.685 sec. 

a. Average Front Wall Overpressure vs Time. Step 1. For P^^ * 
10 psi from figure 3*20, ^pgfl “ ^5*3 psi. 

Step 2. From figure 3.35> h' = 15 ft. 

Step 3* For P^^ » 10 psi from figure 3.21, = 1290 fps. 

Step 4. t^ = 3h'/c^g^^ = 3(15)/1290 = 0.0349 sec. 

Step 5. For P = 10 psi, from figure 3.23> = 2.23. 

SO o 

Steps 6 through 8. In table following. 


From figure 3*12 P^ 


10 psi. 

» 0.262 


Data: I8 KT typical airburst. 

3000 ft distance from ground 


Exterior dimensions of prob- 
lem stmxcture (figure 3 •35). 


3~10b 


t 




0 

0 

1.000 

2.23 

0.0685 

0.10 

0.634 

1.4o 

0.137 

0.20 

0.397 

0.89 

0.2055 

0.30 

0.245 

0.55 

0.274 

0.40 

0.148 

0.33 

0.34-25 

0.50 

0.087 

0.19 

0;4ll 

0.60 

0.049 

0.11 

0.4795 

0.70 

0.026 

0.06 

0..548 

0.80 

0.012 

0.03 

0.6165 

0.90 

0i004 

0.01 

0.685 

1.00 

0.000 

0.00 


p /p 

s' so 
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P P^ + 0.85q 

•S s 

1.000 

10.00 

12.12 

0.8l4 

8.15 

9.48 

0.655 

6.55 

7.40 

0.519 

5.19 

5.71 

0.402 

4.02 

4.33 

0.303 

3.03 

3.21 

0.220 

2.20 

2.30 

0.149 

1.49 

1.55 

0.000 

0.90 

0.93 

0.04l 

0.4l 

0.42 

0.000 

0.00 

0.00 


Saii5)le calcvQatlon of typical entry in table: 


For t/t^ 


0.40(0.685) 


* 0.274. 

From table 3*2 for t/t^ 

= 0.40, <i/q[^ = 0.148 

(j = 0.148(2.23) = 0.33 psi 

From table 3*1 for "^Ao 

= 0 . 40 , P = 0.402 
^ S' so 

= 0.402(10) » 4.02 psl 

P = 0.85q = 4.02 + 0.85(0.33) 
s 

= 4.33 psi 

Step 9. Plot the curve as 
shown in figure 3 *36(81). 

b. Average Back Wall Over- 
pressure vs Time. Step 1 . For 
P = 10 psi, from figure 3.9, U 

SO ^ 

= i 403 fps. 

Step 2 .. 1 ( 3 ^ ■ * 5 Vi ^®3 

= 0.0385 sec. 

Step 3 * \ 

» 4 ( 15)/1115 » 0.538 sec. 

Step 4 . For (t - t(3^)AQ 
= 0.0538/0.685 = 0.0785, from 

table 3.1, PsAso * 

Pg = Pg^ » 10(0.852) » 8.52 psi. 


25e 


^ 20 


o. 



Time (tec) 


(a) Average Overpretture 
vs Time Curve 



(b)Average Net Horizontal 
Overpressure vs Time 
Curve 



(c) Local Roof Overpressure 
vs Time Curve 


(d) Average Roof 

Overpressure vs Time 
Curve 


Figure 3.36. Overpressure vs time curves- 
closed rectangular structure 
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Step 5. 


^^back^max 


p I 1 + ^,l - P^e ^]/2 where P 
0.3^ 

= 8.52 [1 + 0.47] /2 » 6.26 psi 


_ - . [1 + (1 - P)® 

hack‘d max so l. 

p = o.5(io)/i^*7 


O.SP^yi^.T. 


i1 


Steps 6 and 7. 
t - t- 


In table following. 


0.0923 

0.1070 

0.1755 

0.2MK) 

0.3125 

0.3810 

0.1^i^95 

0.5180 

0.5865 

0.6550 

0.7235 


0.0538 

0.0685 

0.137 

0.2055 

0.27^ 

0.3^25 

0.1<-11 

0.4795 

0.548 

0.6165 

0.685 


0.079 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 




so 


0.850 

0.8l4 

0.655 

0.519 

0.402 

0.303 

0.220 

0.149 

0.090 

o.o4l 

0 


8.50 

8.l4 

6.55 

5.19 

4.02 

3.03 

2.20 
1.49 
0.90 
0.4l 

0 


^back'^^s 

0.735 

0.735 

0.740 

0.750 

0.767 

0.790 

0.820 

0.855 

0.897 

0.947 

1.000 


pack 

6.26 
5.98 
4.85 
3.89 
3.08 
2.39 
1.80 

1.27 
0.81 
0.39 

0 


Sample calculation of typical entry in table above: 

For (t - \)Ao = 0 * 5 > 1 = 0.5(0.685) + 0.0385 = 0.3810 
From table 3.1° for (t - t^)/t^ = 0.5, P3/P30 = 0.303 

P„ = 0.303(10) 


"d^' -o 

3.03 psi 


For (t - t^)/t^ =0.5 
^back/^s ~ ^^back^max/^sb ^ ^^back^max/ sb 






t - 


^ J 

= 6.26/8.52 +[i - 6.26/8. 52] [°*^®^o.6B5°-°Co53B°*°^'^'^^ 

= 6.26/8.52 t [1 - 6.26/8.52] [0.2887/0.6312]^ 

= 0.790 

= 3.03(0.79) = 2.39 psi 


back 

step 8. Plot the curve as shown in figure 3.36(a). 
c. Average Ket Horizontal Overpressure vs Time. Step 1. 


Determine 


values of P 


front 


and P 


back 


at times t = 0.00, 0.05, 0.15, 0.20, 


0.65, 0.70, from figure 3.36(a) and determine P^^^ = ^front 


- P, 


"Kq r»"lr 


for* 


each time. 

Step 2. Plot curve of P^^^ vs time as shown in figure 3.36(b). 
d. Local Roof Overpressxrre vs Time at Point x of Figure 3.35« 
Step 1. t^ = L'/Uq = 27/1403 = 0.0192 sec. 
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Step 2. V = (o.042 + 0.108 = [0.042 + O.IO8 (O.5)] l403 

= 134. T fps* 

step 3. t . L'/v = 27/1311.7 = 0.200 sec. 

Step ll. At tlcne t = t_^, = MP^^/Hi.T) ^ 

P JT = Ml0/l‘l.7)(27/5** - 1) + 1 - -0.36 < 0. 

roof' s 

Therefore, = 0 ^ 

p /P = ~ ^ 151 i'/Uq 

steps 5 ttooSi 8. In tetle following. (Note that in the computa- 
tions in table all tim«s in coluhMS (l) and (2) are calculatei from ass^d 
values in column (3) except at times t = t^ - 0.200, t - g (t^ 2 


(0.0192 

t 


0.200) = 0.1096, and t 


0.0192 

O.O8TT 

0.1096 

0.1562 

0.200 

0.22lt.7 

0.2932 

0.3602 

0.3617 

0.4302 

0.4987 

0.5672 

0.6357 

0.7042 


0 

0.0685 

0.0904 

0.137 

0.2055 

0.274 

0.3410 

0.3425 

0.411 

0.4795 

0.548 

0.6165 

0.685 


0 

0.1 

0.132 

0.2 

0.3 

0.4 

0.498 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 


0.3602.) 
s' o 


1.000 

0.814 

0.761 

0.655 

0.519 

0.402 

0.305 

0.303 

0.220 

0.14'9 

0.090 

0.041 

0 


^s 


p 

■^roof 

10.00 

1.0 

10.00 

8.14 

1.0 

8.l4 

7.61 

1.0 

7.61 

6.55 

0.483 

3.16 


0 

0 

5.19 

0.154 

0.80 

4.02 

0.583 

2.34 

3.05 

1.0 

3.05 

3.03 

1.0 

3.03 

2.20 

1.0 

2.20 

1.49 

1.0 

1.49 

0.90 

1.0 

0.90 

0.4l 

1.0 

0.41 

0 

1.0 

0 


Sample calc\alation of typical entry in table above: 

For (t - t^)/t^ = 0.4, t = 0.4(0.685) + O.OI92 = 0.2932 sec 

From table 3-1, ^s^^so " 

P = 0.402(10.0) = 4.02 psi 

Since % J-P varies linearly from 0 at t = 0.200 to 1.0 at 
1*001 s 

t = 0.3602 

p /p = (0.2932 - 0.200)7(0.360 - 0.200) = 0.583 

roof' s , 

P ^ = 0.583(4.02) = 2.34 psi 

step 7. Plot the curve as shown in figure SoolCy- 

e. Average Roof Overpressure vs Time. Step 1. t^ » 

= 54/2(1403) » 0.0194 = 0.0194 sec. 

Step 2. At time t « O.O388 sec 
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P /p . 0.9 + 0. 1(1.0 - p /ii^-7f » 0.9 + 0. 1(1.0 - 10/14.7)^ - 0.910 

roof' s °° \ 

(i.e., this value may not exceed one). 

Step 3. At time t ■ 5 L/Uq » O.1928 sec. Compute 2.0 - 
(1 " fe) f * iBr) - 0-902 and 0.5 . 

0.125 ^2 - 55^^ - 0.5 0.125 ^2 - 5^) - 0.717 

-Since 0.717 < 0.902, "rooA = 

Step 4 . At time t ■ 51 '/Uq t 15^1*/^^ * 0.1928 + 0.l6o4 *• 0.3532 sec 


P Jv = 1. 

roor s 

Steps 5 throvigh 7 - In table following. 


t 

t - t , 
d 

(t - t^)/t^ 

p /p 
s' so 

P 

s 

^roof'^^s 

p 

roof 

0.0388 

0.0194 

0.0283 

0.944 

9.44 

0.910 

8.58 

0.0879 

0.0685 

0.1 

0.814 

8 .l 4 

0.848 

6.90 

0.1176 

0.137 

0-2 

0.655 

6.55 

0.762 

1^.99 

0.1928 

0.1734 

0.252 

0.581 

5.81 

0.717 

4.17 

0.2249 

0.2055 

0.3 

0.519 

5.19 

0.774 

4.01 

0.2934 

0.274 

0.4 

0.402 

4.02 

0.894 

3.60 

0.3619 

0.3425 

0.5 

0.303 

3.03 

1.0 

3.03 

0.4304 

0.411 

0.6 

0.220 

2.20 

1.0 

2.20 

0.4989 

0.4795 

0.7 

0.149 

1.49 

1.0 

1.49 

0.5674 

0.548 

0.8 

0.090 

0.90 

1.0 

0.90 

0.6359 

0.6165 

0.9 

0 . 04 l 

0 . 4 l 

1.0 

0 . 4 l 

0.7044 

0.685 

1.0 

0 

0 

1.0 

0 


Sample calculation of typical entry in table above: For (t - 

= 0.1, t - 0.1(0.685) + 0.0194 - 0.0879 sec 

From table 3-1 for (t - t.)/t « 0 . 1 , P /P « 0 . 8 l 4 
' d'' o ' s' so 

P * 0 . 8 l 4 (l 0 ) = 8.14 psi 

S 

Since Pj.QQf/Pg varies linearly from O.9IO at t ■ O.O388 sec to 
0.717 at t = 0.1928 sec 


’’r-oo/^a " : oioii ) 

= 0.910 - 0.0491(0. 193 )/o. 1540 
= 0.848 


roof 


= 0 . 848 ( 8 . 14 ) = 6.90 psi 


Step 8 . Plot the curve as shown in figure 3 * 38 ( 1 ). 
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TmPHTG ON HECTAlTCiUIAR STRUCTURES WITH OPEMINGg 

3_11 WFEOTS OF OFEMINGS ON LOADING, For the determination of the loads 
on rectangular structures with openings it is necessary to consider all 

structures as falling into one of two categories. 

Category (l). Structures having open interiors free of walls and 
other obstructions so as to allow relatively unhindered propagation of the 
■blast vave through the Interior of the structure, once it has entered 
throu^ windows or otherwise. It is assumed that the percentage of hack 
face openings is approximately equal to the percentage of front face 

openings . 

Category (2). Structures having interior walls and other construc- 
tion which would prevent free propagation of the blast wave through the 
interior of the structure. 

Those structures falling into Category (2), above, are analyzed in 
a manner similar to that for the case of closed rectangular structures when 
determining the net translational forces acting on the structure. Net 
loads on exterior and interior portions of these structures may be computed 
as are the loads for structures falling into Category (l), vith suitable 
inteirretation of the various structural dimensions involved. The deter- 
mination of the loads on structures with relatively open interiors. Cate- 
gory (1), is described below for approximately equal percentages of open- 
ings in the front and. back walls . 

The sequence of phenomena fo3.1oWing the impingement of a shock wave 
on a rectangular structure with openings can be briefly summarized in the 
following paragraphs. 

When the shock front strikes and reflects from the front wall, the 
overpressure on that wall rises ‘to the reflected overpressure. Rarefaction 
waves immediately move in from the outside edges of the wall, clearing the 

reflected overpressvires. 

Windows and doors in the front wall will probably break under the re- 
flected overpressure, thou^ not Instantaneously. However, they gen y 
will break before the clearing is complete; hence, part of the clearing 
will occur throu^ the window openings, allowing high pressure air to flow 
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into the interior. This sudden release of hi^ pressure will cause shocks 
to form in the interior. These inside shocks, formed at each opening, will 
spread downstream from each opening and will tend to combine into a single 
shock, front. This interior shock is, initially, weaker than the incident 
shock. 

Meantime the Incident shock front has moved past the front wall and 
is sweeping across the exterior roof and sides. The shock then moves 
across the exterior of the back wall of the structure from the sides and 
from above, building up the overpressure on that surface. 

For those structures which have no interior partitions to obstruct 
the passage of the shock throu^ the interior, the shocks entering at the 
front openings will move through the structure raising the overpressures on 
the inside surfaces as they are covered. 

a. Average Exterior Front Wall Overpress'ure ^ppont' time re- 

quired for reflected overpresstires to clear on the front of a solid wall 
struck by a shock wave is expressed in multiples of the time necessary foi; 
a rarefaction wave to sweep the wall once, When walls with 

numerous openings are considered, however, clearing can take place around 
the edges of the openings. This necessitates a modification of the method 
of determining clearing times because the rarefaction wave travels a much 
shorter distance to cover the wall. 

The distance h^ is introduced as the weighted average distance the 
rarefaction wave must travel to cover the wall once, assuming immediate 

access of the Incident shock to the interior 
of the structure. As shown in figure 3*37 
the front face is divided into rectangular 
areas, determined by the location and dimen- 
sions of openings considering the directions 
along which the area can be cleared in the 




Figure 3.37. Subdivision of a 
typical wall with openings to 
determine weighted average 
clearing height 


shortest possible time. These areas are 
labeled in the following manner: 

Step 1. Label 1 , all areas cleared 

SL 

from two opposite sides. 


Step 2. label 1^, all areas cleared from two adjacent sides. 
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Step 3. Label 2 , all areas cleared from orte side. 
Step k. Label 3 , all remaining areas. 

The wei^ted average clearing hei^t is then: 


where 




6 h A 


(3.16) 


A = net area of front face (wb less openings) 

K = area of each portion of the subdivided front face, except 
openings 

h = for areas 1 , the average distance between the sides from 
^ which clearing occurs 

= for areas 1 ^ and 3 , the average hei^t or width, whichever is 
smaller 

= for areas 2 , the average distance between the side from which 
clearing occurs and the side opposite 


5 ■ the area clearing factor, 

n 

a 1/2 for areas 1^^ 
a 1 for areas 1^ 
a 1 for areas 2 


a 1 for areas 3 ^n^n^n 

Z = the notation to represent the summation of the j- 
1:61111 for all areas 

h' = h or w/ 2 , whichever is the smaller 
The tJjK reauired to clear the front face of a structure with open- 


ings of reflection effects t^ 


can now be evaluated from the relation 


t' 

c 


3^f/‘^refl 


(3.17) 


Windows in a structure will not he lasediately blown outi CBEESHOUSE 
data (see page 72 , reference 15 ) Indicates that the wlndow-hreahlng time Is 
Of the order of a millisecond or less. This thre, however, is too short to 
he of any practical tnl)ortance In computing the loading on the structure. 
Therefore, it Is assumed that the windows fail as soon as they are stru 
by the shock front and do not affect the loading computation. 

The time curve of the average outside front wall overpressure on the 
net waU area can he constructed hy the methods outlined in paragraph 3-08a 
and illustrated in figure 3-25 using the clearing time as defined by 
equation ( 3 . 1 ?) • 
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b. Average Interior Front Wall Overpressure ^jL_fpQnt' average 

interior overpressure on the inside of the front wall R ^ ^ initially 

i- front 

builds up in a s imi l a r m a nn er as the average overpressures on the exterior 

back wall of a closed structure, except the vortices are located all around 

the inside edges of all openings. The build-up is linear from zero at time 

t » 0 to the peak overpressure ( see paragraph 3-lle and figures 

3.'41 and 3-^2) of the interior shock at time t « 4h!Vc , where h' is 

if' o^ if 

the wel^ted average build-up 

P*® A height of the interior front 

/ \ wall coin: 5 )uted from equation 

I \ (3*l6) but with all quantities 

p interpreted as applying to the 

S «oi i I interior suirface of the front 

t wall. The overpressure is 

° maintained and builds up again 

ZI] J , . to when the interior shock is re- 

4h'if/co L|/ci Jh'f/Co 

— : : — -J , fleeted back from the rear 


^lo ®o wall to a peak value which is 

Figure 3.38, Average interior front wall overpressure dependent on the size of the 
vs time-small rectangular structure (L./U. ) <0,lt 

‘ " structure. Equlllbrivun is 

then attained when the overpressuire within the interior reaches P : P 

s ^ s 

being the blast wave over- 

p 

pressure existing on the *® \ 

outside of the structure. I \ p 

Figure 3-38 Illustrates the \ 

jf I 

overpressure vs time cxirve 3 / ■ 

for the inside front face s- / I 
of those structures in which ° j 

L./U. < 0.1 t , and \ , ^ ~ ^ . 

^ o' W/co T4 /c Ti™.., *< 

figure 3-39 for those in [*— — ^ 

/ I-|/U|q +L.|/ca 

’Which L./u. > 0.1 t . 

^ ^ ^ Figure 3»39» Average interior front wall overpressure vs 

Average Net Front time^^larse rectangular xtrtmturfi ft. /h 0 It 


Figure 3.38, Average interior front wall overpressure 
vs time— small rectangular structure (L-AJ^^)< O.lt^ 


4 h'lf /co 3 h'f /co 
^ ^ 
L|/Uio+L|/Co 


Time, t 


Figure 3,39. Average interior front wall overpressure vs 
time— large rectangular structure (L./lj^^)> 0.1 1^ 


Wall Overpressure average net overpressure acting on the 
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\ ^--Average Exterior Front Woll OverpressurOi Pf^ont 
^-Average Interior Front Wall Overpressure, Pj. 
✓Average Net Front Woll Overpressure, Pf_n,f 


Time, t 

Figure 3 » 40 » Average net front wall overpressure 
vs time— rectangular structures with openings 


front wall is deter- 

mined by subtracting the Interior 

front wall overpressure from the U _ _ 

\ ^ “ ^front “ ^1- front 

exterior front wall overpressure \ \ 

on a common time basis according g ^^^Av«rog« Exterior Front Woll 0v«rpr«»»ur6, 

to the equation | X'j^^Averoge interior Front Woll Overpreesure. 

— — — - Q\ o ✓Average Net Front Woll Overpressure, Pf_ 

^f-net ^front " ^i-front ^3* ; I 

— I — 1 \ ij- 

where P . t P« and / \/v 

wnere front^ / y — 

p are given in terms of 

i-front , 

t. ngure 3.40 illustrates 

graphically the variation of the Figure 3 . 40 . Average net front wall overpressure 

vs time— rectangular structures with openings 
average net front wall overpres- ® 

sure with time. 

d. Average Exterior Back Wall Oveipressure ^i 3 a,ck* ^ "tin's lag is 
present between build-up of pressiire on the back wall due to the time re- 
quired for the blast wave to travel the length of the structure. Until the 
•,jriii(3_ows or window frames in the back wall fail^ the back wall loads behave 
in the same manner as those on the back wall of a closed structure. After 
the windows break the back wall loading is affected by the vortices formed 
at all of the edges of the rear openings due to the interior blast passing 
through the openings at the same time that the exterior wave sweeps the 
rear face. If it is assumed that the windows blow out instantaneously over 
the entire stiucture, l.e., the window-breaking times are zero^ the build- 
up time for pressures over the rear face is 

H = 4h^/c^ (3.19) 

vhere is the weighted average build-up time of the rear face, defined 




5 h A 
n n n < . , 

\ " 


( 3 . 20 ) 


Apj is the net area of the back wall (gross area less openings) 

6 , h , A are as defined in paragraph 3-lla 
n^ n^ n 

h/ is computed in the same maimer as h|. (the weighted average 
^ clearing height of the front face) 

The average overpressure vs time curve for the outside rear face can 
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be coniputed by the methods outlined in paragraph 3-08b and illustrated in 
figure 3.27 using the build-up times as given in equation (3.I9). 

e. Average Interior Back Wall Overpressure P^. The presence 

of openings' in the front vail of a structure pennits^I portion of the inci- 
dent shock to enter through these openings, producing an interior shock 
wave weaker than the incident, which becomes approximately plane as it 
travels the length of the interior. That portion of the interior shock not 
passing through the back wall openings is reflected by the back wall, ac- 
cording to the laws of normal reflection and returns toward the front. 
Q-aick clearing occurs on the inside and vortices form at edges of openings 
on the outside as the blast wave passes out through openings. 

In the prediction of inside pressures the structure is assumed to be 
smooth. Structural features which do not follow this idealization are 
neglected as regards their influence on the determination of Interior 
loads. These features include interior columns along side walls, pipe 
ctions which cross the interior, etc. The peak overpressure of the in- 
shock wave is dependent upon the overpressure of the incident 

shock wave and the percentage of window openings. Figure 3-41 is a plot of 

soi ^SO different values of the ratio of the area of openings in 

the front face to the gross area of the front face Fl^e 3-42 

s a plot of vs for various peak incident overpressures, 

^ese curves are based upon data reduced by the Armour Eesearch Foundation 

L 7] from shock tube tests performed at the University of Michigan [16] and 
at Princeton University [l8] . 

tte inside back remains nnloaded nntil time t =. b /U. , where 
Lj Is the distance from the outside of the front wall to the insSe rear 
face and is the velocity of the Inside shook front which Is plotted 

in figure 3.9. The overpressure on this wall Is then Instantaneously 
raised to the reflected overpressure valrre given m figure 3.a0. 

Ihe presence of openings in the back wall prevents the maintenance of 
this reflected overpressure and, hence, there Is a linear decay to a value 
approachl^ the drag pressure In the Inside blast wave. The time required 

wall 1° nr 

o the sise of the structure. For relatively small structures 
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Peak Incident Overpressure, Pg^lpsi) 

'igure 3.41. Peak interior overpressure vs peak incident overpressure for various percentages of front wall openings. 

(Validity of dashed portions of curves uncertain) 
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there is an additional build-up following an initial decay due to the re- 
turn of the shock wave reflected from the rear wall toward the front wall, 
while for longer structures this second build-up is absent. It is assumed 
that the interior blast wave eventually reaches a condition of equilibrium 
wherein the interior and 


the incident blast waves 
have eq,ual overpressures 
and the back wall inte- 
rior overpressure eq.iials 
the drag pressure in the 
incident blast wave, with 
the time displacement 

factor t^ = ^i^^io* 

To acco\ant for the 

effect of size of the 
structure on interior 
overpressures, two curves 
of the interior back wall 
overpressvires are given. 
Figure 3.^^3 is applicable 
to structures in which the 




Figure 3.44. Average interior back wall overpressure vs 
time— large rectangular structure ^ ^"^*0 

6l 


Figure 3.43. Average interior back wall overpressure vs 
time-small rectangular structure (L/ V.J < O.lt^ 


time required for the inside 
reflected shock wave to 
travel the length of the in- 
terior is less than one- 
tenth the duration of the 
positive phase of the out- 
side shock wave; i.e., 

Li/Uio< 3*^^ 

is applicable to those struc- 
tures for which > 

0 It . P , •, as a func- 

u.xXq* -^i-back 

tion of time is as illus- 
trated in these figures. 
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f . Average Net 
Bask Wall Overpres- 


sure 


b-net* 


The 


average net overpres- 
sure acting on the 

back wall P, . is 
b-net 

determined by sub- 
tracting the exterior 
back wall overpres- 
sure from the inte- 
rior back wall over- 
pressure on a common 
time basis according 
to the equation 

^b-net “ \-back “ ^back (3*2l) 

\ack 4 . Figure 3.45 

illustrates graphically the variation of the average net back wall over- 
pressure with time. 


Figure 3.45. Average net back wall overpressure vs time- 
rectangular structure with openings 


where 


^b-net' ^i-back^ 


g- Average Exterior Roof Oveipressure In this paragraph it 

is asstuned that the roof is plane, horizontal, and without roof surface 
openings. Therefore, the average overpressure vs time cmves on the ex- 
terior are computed by the methods outlined in paragraph 3-08e and illus- 
trated in figure 3.34. The quantity h' used in these figures should be 
replaced by h^ . 

h. Average Interior Roof Overpressure Here again it is 

necessary to place the structure being considered into one of the two cate- 
gories established in the previous paragraphs due to the difference in be- 
havior of the inside shock. For the smaller interiors the average over- 
pressure on the inside of the roof will build up in time L^/u to the 

initial peak overpressure of the Inside shock wave, at which ti^e reflec- 
tion of the interior shock from the rear face will take place. This re- 
flected shock then travels toward the front wall spreading higher pres- 
sures along the underside of the roof. The reflected shock front is 
destroyed by the vortices at the front face as well as by flow out through 
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openings and no fur- 
ther reflections oc- 
cur. The average roof 
overpressiire then de- 
•creases from the rear 
face reflected value 
to the equilihrium 
value of P_, the In- 
cident air blast 
overpre s sure . Figure 
3 . 14-6 displays the 
average interior roof 
overpre ssiire for 
those structTires in 
vhich 



Figure 3.46, Average interior roof overpressure vs time-^small 
rectangular structure (L./ 


< O.lt . 

o 



Figure 3.47. Average interior roof overpressure vs time— large 
rectangular structure (L./ > O.lt^ 


Figure 3*^7 
plays the average in- 
terior roof overpres- 
sure curve for those 
structures in vhich 

these structures the 

average overpressure 

does not build up to 

P_, since the 

i-refl 

shock wa.''e reflected 
from the rear face 
decays before it can 


spread its higher pressures over the whole of the imderside of the roof. 

1. Average Net Roof Overpressure average net over- 

pressure acting down on the roof ^^_net determined by subtracting the 
average interior roof overpressure from the average exterior roof 


63 


EM 1110-3^^5-1^3 
1 July 59 


3-llJ 



^ ^ . * average 

lirterlor sia« »bij. ovBn,resmjre Is bandied In axactOy the etuK as 

the average Interior roof oveijreeeure. It Is assumed that the side is 

vlthout openings, as for the exterioT* «^rlsa v.«t, 

1- Average Net Side Van 


Overpressure P 


’s-net‘ average 


i^et oveirressure acting inward on 
the side wall is deter- 3 

inined by subtracting the average in- I 
terlor side wall overpressure from ° 
the average exterior side wall over- 
pressure on a common time basis ac- 
cording to the equation 

^s-net * ^side " ^i-side (3.23) 
Figure 3.1*-9 illustrate# graphically 

the variation of the average net 

side wall overpressure with time. 
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3-12 EROCBDURE FOR COMFUTATIOIT OF LOADS ON REC!IAKGULAB STRUCTURES WITS 
OPENINGS. Repeat Steps 1, 2, and 3 luider paragraph 3-09. 


a. 

termine 


Average Exterior Front Wall Overpressure vs Time. 

h’ from the relation 
f 




•L 


8 h A 
n n n 


Step 1. 


De- 


Step 2. Determine h', the clearing hel^t for the front face, treat' 
ing it as if it were closed. Use h^ » h', if Step 1 gives h^ > h'. 

Step 3* Determine t^ from the relation t^ = 

Step 4. Use Steps 1, 2, and 6 to 9 of paragraph 3-09a to obtain the 
^front ™ 

b. Average Interior Front Wall Overpressure vs Time, Step 1. De- 
termine ■til® 1 - ghted average build-up hel^t for the inside front 

wall. 


Step 2. Determine knowing and k^^/k^, from figure 


so 


of'“gf^ 


3 . 1 i-l or figure 3-l*-2. 

Step 3. Determine U^^^, knowing from figure 3 . 9 . 

Step 4. Determine knowing Pg^j^y from figure 3.20. 

Step 5* Determine the curve of P using figure 3*22 or table 3*1* 

s 

step 6 . Plot the curve of illustrated in 

figures 3*38 oi" 3* 39* 

c. Average Net Front Wall OveiTjressure vs Time. Step 1. Knowing 

P„ and , for a sequence of times t^, t-, t„, etc., determine 

front back ^ 1 ' 2 ^ 3 ' 


■ f-net 


P„ - P, , for each value of t. 
front back 


Step 2. Plot the curve of illustrated in 

figure 3 . 40 . 

d. Average Exterior Back Wall Overpressture vs Time. Step 1. De- 
termine h^ from the relation 

6 h A 

, , n n n 

h-e— 


Step 2. Determine h', the clearing height for the rear face, treat' 

ing it as if it were closed. Use h^ - h’, if Step 1 gives h^ > h*. 

Step 3 . Determine t^ from the relation t^ » 4h^/c^. 

Step 4. Use Steps 1, 2, and 4 to 8 of paragraph 3”09b to obtain the 

curve of P-_„v time, 
back 
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®* Average Interior Back Wall Overpressiire vs Ttmft. step 1. Deter- 
mine the wei^ted average huild-up height for the inside hack wall. 

Step 2. Determine and as in paragraph 3-12b. 

Step 3* Determine t, = L./U. . 

d i' lo 

Step 4. Determine overpressure at time t - t, = L^/c from the 

d i' o 

relation 

^i-hack * ^soi + " ■^of/^gf^ ^^1-refl “ ^soi^ 

Step 5. Determine the curve of P^ + 0.85q as in Steps 1 and 6 to 8 
of paragraph 3-09a. 

Step 6. Plot the curve of P. , vs time as illustrated in 
figure 3.43 or 3-44. 

_ ^ei^ge Net Back W all Overpressure vs Tlmft. step 1. Knowing 

fi-hack \ack_ ^ sequence of times t^, t^, ty etc., determine 

^h-net “ ^i-hack ‘ ^back each^value of t. 

Step 2. Plot the curve of vs time as Illustrated in 

figure 3.45. 

®* Average Exterior Poof Overpressure vs Time. Step 1. Use Steps 
1 to 8 of paragraph 3-09e to determine the. curve of P^^^^ vs time, but 
replacing the quantity h' appearing in any of these Steps by h^. 

garage Interior Poof Overpressure vs TOmp. step 1. Determine 

^soi' ^io' ^i-refl paragraph 3-12b. 

Step 2. Determine t. » L^/2U. . 

d V io 

Step 3. Determine the curve of P^ using figure 3.22 or table 3.1. 

Step 4. Plot the curve of vs time as illustrated in 

figure 3.46 or 3.47. 

Average Net Roof Overpressure vs Tlmp. step 1. Knowing P 

^i-roof of times t.., t-, t., etc., determine P 

mV -"d n, s. d ^ r-net 

roof ’^i-roof value of t. 

step 2. Plot the curve of P as time as Illustrated in 

figure 3.48. 

J* ^rage Ebcterlor Side Wall Overpressure vs step 1. Use 

the procedure of paragr^h 3-09g to determine the curve of P vs time 
but replacing the quantity h' by h^. 
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k. Average Interior Side Wall Overpressure vs Time. The cxirve of 

^i-slde obtained in the same manner as as presented in para- 

graph 3“12h. 

l . Average Het Side Wall Overpressure vs Time. Step 1. Knowing 


and P 


i-side ^ sequence of times t^, t^, t^, etc., determine 


for each value of t. 


s-net side i-side * 

Step 2. Plot the curve of P + vs time as illustrated in 

S*“I1© u 

figure 3<^9* 


LOADIITG OK STRUCTUBES WITH GABLED ROOFS 

3“13 LOADING ON GABLED ROOFS, a. Average Front Slope Roof Overpressure 

P Shock Front Parallel to Ridge Line. When the shock wave strikes a 

roof 

roof surface, or the earth cover of 
a semiburied structure, pitched as 

in figure 3.50(a) or (b), oblique bim^ l 

reflection of the shock front occurs. LmmiimiA. 

The reflected pressure is a 

function of the peak incident over- ► 

pressure P „ and the angle of in- [T l II l 11 1 

clination of the roof 0. The ratio 


P^ q/^so plotted in figiire 
3.11 as a function of the angle of 



Figure 3.50. Typical gabled roofs 


incidence a, here equal to 90 “ ®* 

For a given angle of inclination of the roof, the ratio of the maxi- 
mum average overpressure to the reflected overpressure is a constant when 
the overpressures are expressed in terms of their deviation from the over- 
pressure in the incident shock wave. This relation is expressed as follows, 
where Kg is independent of the peak incident overpressure 


where 


P’ 

roof 


peak average oveipressxire on front slope of each bay at 
time t « t + lAu 


length of one gabled roof span 


lAu 
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r-a 


= reflected overpressure as obtained from figure 3. 11 



The factor is plotted as 

a function of 0 in figure 3.51. 
The peak average overpressxare on 
the front roof slope P* ^ is 
then computed from the following 
equation 

^roof “ ^sL ^ ^e^^r-a " ^sL^ 

(3.25) 


0 ® 10 ® 20 ® 30 ® 40 ® 50 ® 60 ® 70 ® 

Inclination of Roof, ^ 

Figure 3,5L Ratio defining peak average 
roof overpressure vs inclination of roof 

by equation (3.25) at time 
clears to the 
the relation 


The average overpressure on 
the front slope of the roof builds 
up linearly from zero at time 
t = t^ - to the value given 

The average oveipressure then 
= t^ + 3 L/ 4 u^ and is given by 


t = t , + lAu . 

CL O 

drag overpressure at time t 


for all times 


in excess of t 


P - ® P + C^q 
roof s 

“ t^ + 3 LAu^^ where 


(3.26) 


Proof “ average overpressure on 
the front roof slope 

Pg « overpressure in the inci- 
dent shock wave 

^d “ time-displacement factor 
» pAu for front slope 
of first bay 

q * dynamic pressixre 

Cj^ » drag coefficient plotted 
in figure 3.52 as a 
function of 0 [22]. 

The local overpressure vs time 



Figure 3,52. Drag coefficient for gabled 
roofs vs slope angle of roof 


curve at a point on the front slope of the gabled roof surface consists of 


an instantaneous overpressure rise to the reflected pressure p fol- 
lowed by a linear decrease to the drag overpressure, given by equation 
(3.26), in a clearing time equal to l/ 2 U^. The time displacement for the 
drag pressure curve is equal to the time required for the shock to travel 


from the front edge of the structure to the point at which it is desired 
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to deteiTDlne the local oveiTpressure. 

The local overpressvire vs time curve at a point on the rear roof 
slope consists slnply of the drag overpressure given by equation (3.26). 

The time displacement t^ for the drag pressure curve is equal to the time 
required for the shock to tiavel from the front edge of the structure to 


the point at which it is desired to 
Figures 3.53 and 3*54 
illustrate aveiage overpressure 
vs time curves for the front 
roof slope in the first bay and 1 
In any bay. The first bay time- 
dlsplacement factor t , is | 

given above, and for any other i 

O 

bay it is the time required for 
the shock f jront to travel from 
the front wall of the structure 
to a point halfway across the 
front roof slope of the bay. 


determine the local overpressure. 

Note: td*L/4Uo L Is Minimum Length of Roof 



Time, t 

Figure 3,53. Average overpressure on windward 
roof slope of first bay of gabled structure vs 
time^^ridge line parallel to shock front 


P roof 




a 

I. 

m 

> 

O 



!<!+■ tp 


Figure 3,54. Average overpressure on windward roof 
slope of any bay of a gabled structure vs time- 
ridge line parallel to shock front 


The local overpressure 
vs time curve at a point on 
the front slope of the gabled 
roof surface consists of an 
instantaneous overpressure 
rise to the reflected pres- 
sure P , followed by a 
linear decrease to the drag 
overpressure, given by equa- 
tion (3*26), in a clearing 
time equal to The 

time displacement for the 
drag pressure curve is equal 
to the time required for the 
shock to travel from the 
front edge of the structure 
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1*^ P +00 
. 'to 


+ Coq L Is Mlninom Length of 


Roof Parallel to Shock 
Front for which This 
Curve Is Applicable. 


2Uo Uo 


Time, t 


Figure 3.55. Average overpressure on leeward 
roof slope of first bay of gabled structure vs 
time-bridge line parallel to shock front 


to the point at vhich it is desired to determine the local overpressure. 

h. Average Rear Slope Roof Overpressure Shock Front Parallel 


to Ridge Line. No reflection takes place on the rear roof slope > slope B 
of figure 3 •50* ®ie overpressure builds up linearly from zero at time 

^ t « t^ - "the time at vhich 

NOTE: t«,- 3 L/ 4 Uo shock front reaches the roof 

biostI*" 2 peak^ to the overpressure as given 

by equation (3.26) at time t » t , + 

g u 

p +c a *^0. L I. Mininum L.nflth of L/UU . For the rear slope of the 

^ ® • T'v. ] Poroll.l to Shock ° 

3 I A'Ni Fron* for which Thi. first bay the time-displacement 

• 1/ 1 Curvt I. Applteabl*. , 

0. / I factor is t, » 3L/W • !ntie ayer- 

• / ° 

o /I I age overpressure-time curve on this 

u tL L surface is illustrated in figure 

2u; u; * " 

Tlme.t 3.55. 

Figure 3.55. Average overpressure on leeward overpressure vs 

roof slope of first bay of gabled Structure VS time curve at a point on the rear 

time— ridge line parallel to shock front ^ , 

^ ^ roof slope consists shr^ly of the 

drag overpressure given by equation 
(3*26). The time displacement t^ for the drag pressure curve is equal to 
the time required to travel from the front edge of the structure to the 
point at vhich it is desired to determine the local overpressure. 

c* Average Overpressures, Shock Front Perpendicular to Ridge Line. 
The overpressures on all surfaces^ for the blast wave traveling along the 
ridge line of the structiire^ are obtained by the methods outlined in para- 
graph 3 '" 08 . For closed structures^ the clearing height h* is taken as 
the average height of the front vail, or half the vidth, whichever is the 
smaller; the build-up height of the back vail is taken as the average 
hei^t of that vail or half the vidth, whichever is the smaller. The 
clearing and build-up heists for gabled structures with openings are com- 
puted by the procedures in paragraph 3-11. 

3 - 1 ^ procedure f oe computation of loads on gabled roofs. Repeat Steps 1 , 
2, and 3 under paragraph 3-09 . 

^erage Front Slo pe Roof Overpressure In First Bay vs Time, Shock 
Parallel to Ridge Line. Step 1. Determine U^, knowing P , from 
figure 3.9. 
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Step 2. 

Determine 

ta ■ 

Step 3- 

Determine 

Kg, knowing 0, from figure 3. 51. 

Step 4. 

Determine 

PgL^ the overpressure in the incident shock wave 

at time t - t 

tl 

Use figure 3. 10 to get t and figure 3.22 to 


Step 5* 

Determine 

P , knowing a = 90 - 0 and P , from 

r-u so' 

figure 3*11. 

Step 6. 

Determine 

P^oof^ the average maximum overpressure on the 


roof at time t = t , + L/4 u , from the relation 

a o 

P' ^ = K^(P - P ^ ) + P ^ 
roof 0' r-Q! sL sL 

Step 7* Determine knowing 9 , from figure 3-52. 

Step 8. Determine P + by use of Steps 1, 6 , and 7 of para- 

S JJ 

graph 3-090' and the following table. 


time, t t - t^ 

(t - 


q 

Si V^sc S 

(1) (2) 

(3) 


(5) 

00 

vo 

Colvrain (1) 

is a sequence 

of times 

from 

t = t- to t = t. 
d c 


(9) 


Column (4) is obtained from figure 3*24 or table 3 *2 for the values 
in colimin (3)- 

Column (7) is obtained from figure 3*22 or table 3'1 for the values 
in column (3). 

Step 9. Plot the curve of 'vs time as illustrated in 

figure 3'53* 

b. Average Rear Slope Roof Overpressure in First Bay vs Time, Shock 
Front Parallel to Ridge Line. Step 1. Determine U^, knowing P„^, from 

■“ " ' ' -II Q gQ 

figure 3.9. 



Step 

2. 


Step 

3 . 

3-l4a 

» 



Step 

4. 

3.55- 




-r 

s D 


roof 


vs time as illustrated in figiire 
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LOADING ON EXPOSED STRUCTURAL MEMBERS 


3-15 


3-15 I^ADING ON EXPOSED STRUCTURAL mAMENORKB. ae llast wave loading on 
structures with, open steel frameworks such as bridge trusses, steel towers, 
or steel Industrial buildings following the destruction of frangible walls 
is caused almost entirely by the drag pressure resulting from the hi^- 
velocity blast wind which follows directly behind the shock front. The 
load due to the unbalanced pressures resulting from reflection and diffrac- 
tion of the shock front can be neglected because of the extremely short 
duration of these effects. A structure of this type consists of structural 
elements whose dimensions are so small that their front face clearing times 
are less than 10 milliseconds. In addition the menibers are so small that 
they are enveloped by the blast wave in a similar short length of time. 

The total load on the struct-ure due to the drag forces can be obtained by 
con^juting the load on all the individixal elements of the buildings, front 


DRAG COEFFICIENTS, Cq.FOR VARIOUS 
STRUCTURAL CROSS-SECTIONS 


Shapt 

Direction of Wind 
Co 

Shape 

Cq 

□ 

2.0 

P” 

2.0 

1 

I ‘ 

2.0 

0 

i.O 

2.0 

1 1 ^ 

2.0 

t-H* 

i.6 

t-.|r 

2.0 

L 

2.0 


2.2 

n 

1.8 




^ Standard or Wldo FlonQo Soctions 
t Built up Sections, either Riveted or Welded 


DRAG COEFFICIENTS, (Co)„ FOR WINDWARD (UNSHIELDEC 
TRUSS OF BRIDGES 

(Co), is given as o function of the solidity ratio, G<A'/A, 
where A is the totol area included within the limiting 
boundories of the truss. A' is the octuol orea of the 
members in o plane normal to the wind direction. 

For single span trussesi 

6 <0.25 (Co), » 1.8 

G>0.25 (Co),* 1.6 

For multiple spon or very long span trusses.' 

0.00<G<0.20 (Co), *2.0 

0.20<G<0.30 (Co).* 1.8 

0.30<G<0.90 (Co),* 1.6 

0.90<G<I.00 (Co), *2.0 

Figure 3.56. Drag coefficients for various 
structural shapes and trusses 


and back walls and intermediate ele- 
ments, and adding them together on a 
common time basis. The total drag 
overpressure on an element can be 
con^tuted from the formula 

Drag overpressure * Cj^q (3*27) 

where 

Cjj *5 coefficient of drag ob- 
tained from figures 3.56 
and 3 •57* This coeffi- 
cient is applied to the 
individual structural 
elements 

q » dynamic pressure due to 
the air velocity in the 
incident blast wave at any 
time » pu^/2, and is de- 
termined from figures 3.23 
and 3.24 

where L is the dis- 
tance from the front por- 
tion of the structure to 
the element for which the 
load is being computed 
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Figure 3.57. Drag coefficients for leeward trusses 
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Figure 3. 58. Drag overpressure vs 
time on exposed structural element 


Figure 3*58 illustrates the drag 
overpressure-time cxxrve for an exposed 
element located at distance L from 
the front portion of a structure. Hie 
total load on the element at any time 
is obtained by multiplying the drag 
overpressure (at any time) by the pro- 
jected area of the member transverse to 
the direction of travel of the blast 
wave. The total load on the structure 
at any time is obtained by summing the loads on all the individual members, 
based on time t, referenced to the instant the frontmost elements of the 
structure are struck by the shock wave. 

The overpressure-time curve 
for the front wall supporting ele- 
ments of a structure with frangi- 
ble wall covering is shown in 
figure 3 •59* The overpressure on 
the front wall prior to the time 
of failure is determined as pre- 
viously described for the front 
wall of a closed rectangular 
structure. The time of failure is actually a very short time for frangible 
elements co]i5)ared to t^, and the load transmitted to the supporting ele- 
ments prior to failure may be neglected. After the front wall fails the 
stQiportlng elements are subjected to the drag overpressures given by eq^ua- 
tlon (3.27). The overpresstore-time curve for the back wall supporting 
elements is similar to that for the front wall except all times are dis- 



Time, t 




Figure 3,59. Drag overpressure vs time on front 
wall elements supporting frangible wall panels 


placed by the time -displacement factor t^ * where L is the dis- 

tance between the front and back wall. After failure of the back wall, the 
back wall supporting elements are subjected to the drag overpressxires 
given by equation (3^*27). 

EROCEDURE FOR COMFUTATIOH’ OF LOADS OH EXPOSED STRUCTURAL FRAMEWORKS. 
Repeat Steps 1, 2, and 3 of paragraph 3-09 • 
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Step 4. Determine U . knoving P , from figure 3.9. 

o so 

step 5* Determine t^ » time for shock vave to travel from reference 
point to member imder consideration. 

Step 6. Determine from figure 3.56 or figure 3-57. 

Step 7* Determine q. . knowing P , from figure 3*23. 

o so 

step 8. Determine q., knowing for a sequence of times from 

figure 3-24 or table 3 •2. 

Step 9* Evaluate C^q for these times. 

Step 10. Plot the curve of net drag overpressure vs time as illus- 
trated in figure 3*58» 

LOADIFO ON CYLINDRICAL SURFACES 


3-17 LOADING ON CYLINDRICAL ARCH SURFACES, a. Cylindrical Arch Qyerpres- 
sure P Axis Parallel to Shock Front. 

Figure 3*o0a is a definition sketch for 
the notation used in the determination of 
the loads on a cylindrical arch surface. 

Structvires in this category are lllus- 
ti^ited in figure 3-60b. These consist 



X » -^{cos B'-couB) 





Figure 3.60b. Illustration of structures 
with cylindrical arch surfaces 


Figure 3.60a. Definition sketch for 
cylindrical arch notation 


of esqposed arch structures and sur- 
face struct\ires of various shapes 
covered by an earth fill. If the 
arch or earth-fill surface is not 
truly circular, it is approximated 
by an arc of a circle as illustrated 
by dashed lines in figure 3*60b(c), 
and the loads are confuted for the 
equivalent circular surface. Note 
that in some cases a better approxi- 
mation of the ground s-urfade is ob- 
tained with a gabled roof shape, in 
which case the loads are computed as 
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described in paragraph 3-13 • Procedures given in this paragraph determine 
the air blast loading as a ftmction of time on the cylindrical arch surface 
exposed to the air blast, whether it be the sxirface of the axch on the ex- 
posed structure, or the s\irface of the rounded fill over the covered struc- 
ture. With this procedure, zero time is the instant at which the incident 
shock front first strikes the cylindrical arch surface at its intersection 
with the horizontal ground surface, along the line defined by 9 = 6'. 

Figures 3 *61 and 3*^2 illustrate the variation of the local overpres- 
sure time on the surface of a cylindrical arch. Each of these 



5 



Figure 3.61. Local overpressure on front sur~ 
face of cylindrical arch vs time-axis par- 
allel to shock front, 0° $ d'< 0 < 90° 


Tim®, t 

Figure 3.62. Local overpressure on leeward 
surface of cylindrical arch time^^axis 
parallel to shock front, 90° S $ < 180° 

figures is valid only within a certain range of 6 values, which are de- 
termined by the clearing time peculiar to each. These times and the se- 
quence of overpressures within them are defined as follows: 

For 0 < 0* < 0 < 90 (figure 3-6l): The overpressinre at a given 


point rises instantaneously from zero to P at time t = x/u . The 

r-*c% d " o 

^T-a given in figure 3*11 as a function of a, the angle 


value of 

of Incidence of the shock wave. For the cylinder of figure 3.60a, at = 9. 
Clearing of reflection effects occurs at time 


clearing time t^ is given by 


t = t , + t where the 
d c 


The local overpressure at any point after the clearing of reflection 

effects is given by equation (3.29). P^^^ acts normal to the surface of 

the arch. 

76 


(3.29) 


Dynamic Pressure Coefficient, Cq 


3-17a 


EM 1110-345-413 
1 July 59 


vixere q. is the dynamic pressure incorporating the time-displacement 

factor t, 
d 

Cp is the local drag coefficient, -which, is a function of 0 and is 
obtained from figures 3.63a, b, c, and d and 3.64 as discussed 
below 

Pg is the overpressure in the incident blast wave at any time 
t - 



9 (degrees) 


Figure 3.63a. Local dynamic pressure coefficient for cylindrical arch- 
axis parallel to shock front-supersonic Mach numbers (My 1) 
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Figure 3.63h, Local dynamic pressure coefficient for cylindrical arch— axis parallel 
to shock front-^high subsonic Mach numbers (0,4 <M < 1) 


The distribution of dynaioic overpressures about an arch is a fimctioii’ 
of both the Reynolds number and the Mach nimiber M of the high- 

velocity wind in the blast wave [48, 54 J. The Reynolds number, a dimen- 
sionless quantity of the type is- a parameter which characterizes the 

M* 

relative ln^rtance of viscous action in steady nonmlform flow. Since the 
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Figure 3.63d, Local dynamic pressure coefficient for infinitely long cylindrical arch — cucis 
parallel to shock front Q? < 5(10)^ and M < 0.4] 
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Figure 3.64. Mach number of undisturbed stream vs time ratio for various peak overpressures 


Viscosity n and the density p appear as a ratio in the Reynolds number 
it is convenient to treat this ratio of fluid properties as a property 
itself. Thus, the Reynolds number involves only the arch diameter R, a 
velocity u, and the fluid property p/ji. The Mach number, also 
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dimensionless quantity, Is the ratio of the actual velocity to the velocity 
of sound. 

In the following discussion, the distribution of dynamic pressures 
about the circular arch is assumed to be the same as the distribution about 
a circular cylinder having the same diameter. Then the results of e3q)eri- 
mental studies conducted on flow past circular cylinders can be utilized 
for this problem. This procedure will introduce some error but it is be- 
lieved to be insignificant conpared to the mcertainty associated with the 
basic data. 

For flow at low Mach numbers, the pressvire distribution and the total 
drag force change abruptly for a Reynolds nxmiber range between 3(10)^ and 
5 ( 10 )^. This transition is very sensitive to the smoothness of incident 
flow, temperature gradients, and roughness of the surfaces. The explana- 
tion advanced for this change in pressure distribution and decrease in drag 
is that the boundary layer changes from laTn-friftT to turbulent at this crit- 
ical Reynolds number and causes the separation point to move farther back, 
resulting in the redistribution of dynamic pressure observed. 

As the Mach number pf the flow is increased, this decrease in total 
drag and redistribution of pressure at the critical Reynolds n\maiber becomes 
less pronounced \mtil at a Mach number of about O.ll-, it no longer occurs. 
The separation point does not move when the boundary layer changes from 
laminar to t'urbulent and is probably fixed by the presence of shock waves. 
For flows with a free stream Mach number greater than the pressure 

distribution is not influenced by the Reynolds number but it is somewhat 
dependent upon the Mach numiber, the distribution being different for super- 
sonic flow than for subsonic flow. 

Four press\ire distributions are necessary in order that the entire 
ranges of Reynolds and Mach numbers can be considered. The critical values 
of Reynolds and Mach numbers are arbitrarily taken as R * 5(10)^> M = 0.4, 
and 1.0. These dynamic pressure coefficients are applicable for steady 
flow conditions, and since the flow about the arch is far from steady, 
their use may Introduce some additional error. The order of magnitude of 
the error could be ascertained by comparing the results obtained by this 
procedure with direct measurements in field tests in which unsteady effects 
are present. 
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The expression for the local pressure coefficient when R < 5(l0)^ is 

Cjj = ^ [o.T + 0.5 Vi - V(Vd)^ ] (3.30) 

•where 

r' is the local pressure coefficient for an arch infinitely long 
D c r \ 

when R < 5(l0)^ (figure 3-o3<i) 

L is the length of the axis of the arch 

D is the- diameter of the arch 

The Mach number is plotted in figure and the ratio of Reynolds 

number to cylinder diameter is plotted in figure 3-65 for a range of inci- 

dent peak overpressures as a function of ^ 



Figure 3.65. Reynolds number ratio of undisturbed stream vs time 
ratio for various peak overpressures 


Eor 90° < 0 < 180° (figure 3-62); The overpressure at any point on 
the leeward side of the cylinder rises instantaneously from zero to a 
finite overpressure at time t - t^ = 0, given by 

This initial overpressure clears to the value of given hy 

ecLuation (3-29) at time t = t^ + t^, where 

*0 = (V=,efl)(8/9°) 
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Appendix A, Transonic Drag Pressures, appearing at the end of this 
has been extracted from a report prepared by Amman and Whitney, 
Consulting Engineers, under contract DA 49-129-Eng-319 with the Chief of 
Engineers. 

A method is outlined and charts included for evaluating pressures on 
earth- covered aboveground structures that have sloped embankments and a 
limited flat earth fill over their top. The method is based on aerodynamic 
principles which have a laboratory background. The results obtained are 
for a flat- top shock which must be converted to a shock that varies with 
time when applied to the analysis of structures to resist atomic weapons. 

The application to shock waves which vary with time can be accomplished in 
a similar manner described above for the equivalent circle method. At the 
present time the method outlined in appendix A has not been sufficiently 
developed to be set forth as a preferred method of analysis. Supplemental 
data will be available at a future time. 

b. Average End Wall Overpressure Axis Parallel to Shock 

Front. The average overpressure 

P , on the ends of a cylinder 
end 

oriented as in figure 3*60a is given 
in figure 3.66. It is assumed that 
the variation of incident overpres- 
sure in the time required for the 
shock wave to travel across the cyl- 
inder is linear and that the average 
overpressure on the ends is equal to 
the overpressure in the incident 
blast wave P^ with the time-displacement factor t^ = D cos 0'/2U^. 

c. Local Cylindrical Arch Overpressure, Axis Perpendicular to Shock 
Front. The local overpressure at any point on the periphery of the cylin- 
drical arch oriented with its axis perpendicular to the shock front is the 
overpressure in the incident blast wave P^ with a time- displacement fac- 
tor t^ = L’/U^ where L' is the distance from the front end to the point 
at which the overpressure is desired. 



ends of cylindrical arch vs time-axis 
parallel to shock front 
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d. Average End Wall Overpressure, Axis Perpendicul ar to Shock Fron;U 
The average overpressure on the front and hack ends of a cylindrical arch 
oriented with its axis perpendicular to the shock front can he computed as 
outlined in paragraphs 3-08a and 3-08h. The clearing height for the front 
wall and the huild-up height for the hack wall are equal to h as shown 

in figure 3-80a. 

3.18 PBOCEIDUBE for computation of loads on CYLIMDRICAL ARCBES. Repeat 

Steps 1, 2, and 3 of paragraph 3-09* 

a. Local Cylindrical Arch Overpressure vs Time , Axis Parallel, to 

Shock Front. Step 1. Detemine knowing a = 9 and from ^ 

figure 3.11 for 9 < 90°. Initial = (1.5 - e/l80)P^^ for 6 > 90 . 

Step 2. Detemine knowing P^^, from figure 3.21. 

Step 3. Detemine = 4h/c^^^3_ for 9 < 90° or 
(0/90) for 9 >90°. 

Step 4. Determine knowing P^^, from figure 3.9. 

Step 5. Detemine t^ = x/U^. 

Step 6. Detemine at which M = 1 from figure 3-64. 

Step 7. Detemine for times less than this value from figure 

3.63a. t - t , 

Step 8. Detemine at which M = 0.4 from figure 3.64. 

Step 9. Determine for times less than this value hut greater 

than the time detemined in Step 6 from figure 3.63h. 

step 10. Detemine at vhich K = 5(10)^ from figure 3.65. 

o 

step 11. Detemine Cj, for times less than this value but greater 

than the time detemined in Step 8 from figure 3.630. 

Step 12. Detemine for times in excess of the value detemined 

in Step 10 from equation (3.30) and figure 3.63<i. 

Step 13. Detemine P^ + C^q hy following the procedure of Step 8, 

paragraph 3-l4a. 

Step 14. Plot the curve of P^^^ vs time as illustrated in figures 

3.61 and 3.62. 
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Oyerpressure vs Time, Axis Parallel to Shock 
Front . Step 1. Determine t , = D cos 9 ' /2U . 

CL ^ O 

step 2. Determine values of for a sequence of times ranging 

from t - t^ = D cos 0 '/ 2 U^ to t - t^ = t^ from table 3-1 or figure 3 . 22 . 
Step 3. Plot the curves of P^^^ vs time as illustrated in figure 

3.66. 


4 ocal Cylindrica l Arch Overpressure vs Time, Axis Perpendicular 
— Front. Step 1. Determine t^, the time required for the shock 

front to travel from the front end of the arch to the point under con- 
sideration. 

Step 2. Determine P^ for a sequence of times, with a time- 

displacement factor t . 

d 

Step 3. Plot the curve of P = p 

cyl s 

^e?age End Wall Over pressure vs Time, Axis Perpendicular to 

Sh ock Front. Step 1. Determine the clearing height .h' = h (see figure 
3 •63a) . 


Step 2 . Plot the curves of 
methods of paragraphs 3-09a and b. 


^front 

respectively. 


'back 


vs time using 


LOADING ON DOMES 


jading on SPHERICAL DOME SUEFACES. The notation convention used to 
designate a point on the surface of a dome is illustrated in figure 3.67. 
With the notation adopted, any point on the periphery of a spherical dome 
is located by the two angles e and 0 . Denoting by p any point on the 
sphere, 0 is the angle between the horizontal diameter parallel to the 
direction of travel of the shock and the radius Cp joining the surface 
point to the geometrical center of the sphere. The elevation of the point 
P^ above a horizontal plane through the center of the sphere is R sin 0 
sin 52(, where ^ is the angle between the horizontal plane through the 
center of the sphere and the inclined plane containing points 0 , O', and p 
The variation with time of the local oveipressure P, noimal to 
the surface of the dome is similar to that for a cylinder aHhown in 
figures 3.61 and 3.62. The overpressure at a given point rises 
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instantaneously from zero to 
at time t^ = where x = (d/ 2 ) 

- E cos 6 with d being the dis- 
tance across the dcane at its base and 

U is the shock front velocity. The 
o 

angle of incidence at the in^jinging 
shock wave is equal to 9. For 
values of 6 less than 90 ° 'the re- 
flected overpressure is obtained from, 
figure 3 . 11 j for 9 greater than 90 ° 
the wnylmuTn overpressure is given by 
the relation 


'^doma - (^-5 - ®A80)P^^ 


(3-33) 

The initial maximum overpres- 


sures are cleared in time t as 

c 

given by equation (3*34), where h 
is the height of the dome above its 
base. 

= 3 h/c_^, for e < 90° (3.34) 



'refl 

'refl" 


Figure 3.67. Definition sketch for 
spherical dome notation 


t^ = (3h/c_^,)(0/9O) for 0 >90° 


The local overpressure at any point after the clearing of reflected 
overpressures is given by equation ( 3 . 35 ). 


'dome 




(3.35) 


where 


q = unit drag pressure incorporating the time-displacement factor t^ 

= local drag coefficient plotted as a function of 0 in figures 
^ 3.68 and 3.69 


The ratio of the Reynolds number of the air blast wind to the diam- 
eter of the sphere is plotted in figure 3*65 for various incident oveipres- 

tr 

sures as a function of (t - At the Reynolds number E^ = 10^, the 

value of the drag coefficient changes abruptly. The ratio of (t - ■t(^)/t^ 
associated with this value of R^ determines the applicability of either 
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figure 3*68 or 3-69 to oltaln Cj^ for use in equation (3.35). 

The loading in domes calculated by the methods presented in this 
paragraph can be in error for large oveipressures . Use of the present 
method should not be made for pressures greater than I5 psi. 

3-20 PROCEDURE FOR COMPUTATION OF LOADS PIT SPHERICAL DOMES. Repeat Steps 
2, and 3 of paragraph 3-09 . 

Step 4. Determine P^_^, knowing a = 0 and P^^, from figure 3. 11 
for 9 < 90°. Determine maximum P^^^ = (1.5 - 0/l8o)Pg^ for 0 ^90°. 
Step 5. Determine knowing P^^, from figure 3. 21. 


Step 6. Determine 
(0/90) for 0 > 90°. 




for 6 < 90 or 


t = (3h/c 
0 ' refl"^ 


Step 7. Determine U^, knowing P^^, from figure 3.9. 

Step 8. Determine t^ = x/U . 

step 9. Determine (t - t^)/t^ at which R^/d = 10^/D, knowing 
^so^ figure 3.65. 

Step 10. Determine Cj^ for times less than the value of t - t 

obtained in Step 9 from figure 3^68; for times greater than this value 
from figure 3.69. 

Step 11. Determine P^ + C^q by following the procedure of Step 8, 
paragraph 3-l4a. 

Step 12. Plot the curve of P^^^ vs time for the point selected, 
as illustrated in figures 3.61 and 3.62. 


LOADING OH BURIED STRUCTURES 

LOADIIJG OU buried structures. The method for deteimination of earth 
pressure loading on buried structures is more approximate than that for the 
loading on surface structures. The soil has an appreciable mass which af- 
fects the magnitude of the oveipressures acting on a buried structure be- 
cause of the motion of the structural elements which are in contact with 
the soil. This is handled by an empirical procedure in which the load is 
computed independently of the motion of the structural element. However, 
the mass of the element in contact with the soil is increased by the amount 
of the mass of the soil adjacent to the structure, but of a thickness not 
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Figure 3.68. Local drag coefficients for spherical dome, R >10 
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Figure 3.69. Local drag coefficients for spherical domCf R < 10 
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to exceed the span of the element. Other factors which complicate the 
situation are "arching action" and wide variations in the density, moisture 
content, and seismic velocity of the soil with depth and location which 
affect the transmission of pressure waves throu^ the ground. The treat- 
ment below, which neglects "dynamic arching," is largely substantiated by 
the conclusions of reference [46]. 

An air blast wave impinging upon the surface of the groiuad induces 
underground overpressures of a magnitude approximately equal to the blast 
wave oveipressure at shallow depths of 5 to 20 ft [30] in the ground below 
which cause earth oveipressxires on the sides, roof, and bottom of buried 
structures. The pressure wave thus induced in the ground is propagated at 
the velocity C^, the seismic wave velocity in ground (see table 3.3). If 
the ground surface above the buried structures is located in the region of 


Table 3.3. Soil Factors 


Soil Type 

Soil 

Weight 

(Ib/ft^) 

Cg, Seismic 

Velocity 

(fps) 

k. Soil Rres- 
sure Factor 
(Eq (3A0)) 



Min 

Max 

Min 

Max 

Top soil (light dry) 

85 

600 

900 

262 

590 

Top soil (moist, loamy silt) 

95 


1,300 

812 

1,370 

Top soil (clayey) 

Top soil (semiconsolidated 

100 

1,300 

2,000 

1,420 

3,370 

sandy clay) 

110 

1,250 

2,150 

1,510 

4,150 

Wet loam 

Clay (dense wet, depending 

105 

- 

2,500 


5,600 

on depth) 

115 

3,000 

5,900 

8,850 

34,100 

Rubble or gravel 

- 

1,970 

2,600 

6,4oo 

11,100 

Cemented sand 

145 

2,800 

3,200 

9,700 

12,600 

Water-sattirated sand 

125 

- 

4,600 


22, 500 

Sand 

145 

4,600 

8,400 

26,200 

87,000 

Sand clay 

no 

3,200 

3,800 

10,000 

13,900 

Cemented sand clay 

145 

3,800 

4,200 

17,800 

21,700 

Clay, clayey sandstone 

150 


5,900 


45,000 

Loose rock talus 

130 

1,250 

2,500 

1,750 

7,000 

Weather-fractured rock 

160 

1 , 500 

10,000 

3,100 

l 40,000 

Weather -fractured shale 

150 

7,000 

11,000 

63, 000 

156,000 

Weather -fractured sandstone 

160 

4,250 

9,000 

23,500 

116,000 

Granite (slightly seamed) 

170 


10,500 


160,000 

Limestone (massive) 

170 

16, 400 

20,200 

390,000 

590,000 
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regular reflection, the air blast wave will intersect the ground at an 
acute angle, which is assumed equal to zero for computational puiposes. 

The induced ground pressure wave front is then horizontal, nnri moves down- 
ward through the soil and past the structure with a vertical velocity C 
On the other hand, if the ground surface above the stiuctixre is lo- 
cated in the region of Mach reflection, the air shock wave travels parallel 
to the surface of the ground, and the induced ground pressure wave is 
assumed to behave as illustrated in 
figure 3* 70* The leading portion of 
the ground pressure wave is inclined 
downward at the angle y, whose sine 
the surface of the 
ground from the intersection of the 
air shock front with the ground sur- 
face. Although in many cases C 

s 

exceed U^, it is recommended Figure 3.70. Ground pressure wave induced by 
for computation that the angle of wave in region of Mach reflection 

inclination of the ground pressure wave with the earth's surface not exceed 

90° which corresponds to the case where 0^ = 0 ^. As the ground pressure 

wave Impinges on a buried structure, it first loads the upper front comer 
and proceeds horizontally across the roof with a velocity U in the di- 
rection of travel of the air shock wave and simultaneously proceeds verti- 
cally down the front side of the structure with a velocity U x tan y. 

After the ground pressure wave has traveled the length of the^roof, it 
proceeds down the rear side again at velocity U x tan y. 

When the buried structure is located in either the region of regular 
or Mach reflection, the local overpressure on the roof is approximately 
equal to the overpressure of the air blast wave on the earth's surface P 
since apparently no appreciable reflection occurs on the surfaces of a 
buried structure. The lateral overpressure on the vertical sides of a 
buried structure may be considerably less than the pressures applied at 
the top surface of the soil depending upon the type of soil and the height 
of the water table. A value of 0.15P^ has been measured in tests at Nevada 
[46] where the soil was a dry siltv clay. Static tests on soils have 
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demonstrated that the lateral at- rest earth pressure varies from 0.4 to O.5 

of the vertical earth pressure in sandy or granular soils a nd may become 

as large as the vertical earth pressure for soft clay soils because of a 

phenomena called plastic flow [43, 44, and 45]. For the type of dynamic 

loading contemplated here, dry clay soils will probably behave similar to 

granular soils because the load will not persist for a sufficient length of 

time to permit the occurrence of plastic flow. However, it is believed 

that wet clay soils would behave like a viscous fluid and transmit lateral 

pressures instantaneously. For design purposes, it is recommended that 

the lateral overpressure on vertical surfaces of buried structures be 

taken as KP , where K can be evaluated as follows: 
s 

K = 0.25 for dry cohesionless soils 
K = 0*50 for medium cohesive soils 
K = 0.75 for soft cohesive soils 
K = 1.00 below water table 

a. Local Roof Overpressure Regular Reflection Region: The 

local roof overpressure or overpressure at a given point on the roof is 

equal to the air blast wave overpressure on the ground surface P with 

s 

t^ = 0, and where time t is zero when the ground pressure wave strikes 
the roof of the structure. 


Mach Eeflection Region: The local roof overpressure is equal to the 
air blast wave oveipressure on the ground siirface with t^ = 

where time t is zero at the instant the ground pressure wave strikes the 


upper front comer of the structure. 

L * is the distance from the point on 

the roof wheie the pressure is wanted 

to the upper front comer of the 

structure, and U is the air shock 
o 

front velocity. The overpressure- 
time curve is illustrated in 
figure 3.71. 

Average Roof Overpressure 
^roof * I^egu.la-r Reflection Region: 
Average roof oveipressure is 



Tfme, t 

Figure 3.7L Local roof overpressure vs time 
for buried structure located in regions of 
regular and Mach reflection 
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td -L/auo 



td L/U# 


Time, t 


identical with local roof oveipres- 
sure. 

Mach Reflection Region: The 
average roof oveipressure is illus- 
trated in figure 3*72. The over- 
pressure rises linearly from zero 

at time t = 0 to P . the local 

s 

oveipressure existing at the center 
of the roof, at time t = L/U , 
where t^ = L/2U^ and L = length 
of roof in direction of travel of 


Time, t 


r- 9 70 / r . oi zne rooly at time t = L/U 

rigure 3,72, Average roof overpressure vs time ' 

for buried structure located in region of where t^ = L/2U^ and L = length 

Mach reflection 

of roof In direction of travel of 

air shock wave. 

c. Local Front and Back Wall 

Overpressures and P, . . - 

front hack «- S p 

Regular Reflection Region: The local | ‘ 

front and hack wall overpressures are ^ ^ kp, 

equal to KP^. The time-displacement I | 
factor is t^ = h'/u^ tan 7 where 

h’ is the vertical distance from the ° fo+t. 

Time, t 

roof to the point under consideration, figure 3.73. Local front wall, back wall, and 
The overpressure- time curve is illus- side wall overpressure vs time for buried 

tinted in figure 3.73* "r" "T''" 

^ ^ and Mach reflection 

Mach Reflection Region: The local front and rear wall overpressures 

are equal to KP^. The time-displacement factor is t^ = h'/u tan 7 for 

the front face and t^ = L/u^ + h'/u^ tan 7 for the rear face where h' 

^ defined above. The 

I p oveipressure-time curve is illustrated 

” T>C “ 

= / 1 ^ KP, . . 

g KP,o — 7^^'^ -^^srage Front Wall Overpres- 

I ^front’ ^®g^lar Reflection Re- 

° ^ td average front wall over- 

Fi^ur.^ 7 A r pressure is illustrated in figure 3-74. 

tigure3,74. Average front wall overpressure ^p, 

^s time for buried structure located in regions oveipressure rises linearly from 


NOTE: 


2Uotan y 


soh-->^ 


Time, t 


^sgular and Mach reflection 


zero at time t = 0 to KP 
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t = h/U^ tan 7 , where ^ ^ height of the front 

wall of the structure. 

Mach Reflection Region: The average front wall overpressure is illus- 
trated in figure 3.74. The overpressure rises linearly from zero at time 

t = 0 to KP at t = h/U tan y, where t, = h/2U tan y. 
s ' o a o 

e. Average Back Wall Overpres- 


sure 


back* 


Regular Reflection Re- 


gion: The avemge hack wall over- 
pressure is identical with the aver- 
age front face overpressure. 

Mach Reflection Region: The 
average hack wall overpressure is 
illustrated in figure 3 ♦75* The 
overpressure rises linearly from 

zero at time t = L/U to KP 

' o s 

at time t = l/U^ + ti/U^ tan y. 



Time , t 

Figure 3,75» Average back wall overpressure 
vs time for buried structure located in region 
of Mach reflection 


P Regular Reflection Region: 

side 

The local side wall overpressure is identical with the local front and rear 


f. Local Side Wall Overpressure ^ 


wall overpressures. 

Mach Reflection Region; The local side wall overpressure is obtained 
in a similar manner as that for the front and rear wall local overpressures. 
The time-displacement factor is t^ = L'/U^ + h'/U^ tan y where L' is 
the horizontal distance from the upper front comer of the side wall and 
h' is the vertical distance from the upper -front comer of the side wall 


to the point at which it is desired to obtain the overpressure. The local 
side wall overpressure vs time curve is illustrated in figure 3*73* 

g. Average Side Wall Overpressure Regular Reflection Re- 
gion: The average side wall overpressure is Identical with the average 

front face overpressure. 

Mach Reflection Region: The average side wall overpressure is illus- 
trated in figure 3 . 76 . The overpressure rises linearly from zero at time 
t = 0 to KPg at time t = L/U^ + h/U^ tan y, where t^ = (l/U^ + h/U^ 
tan y)/2. 

h. Overpressure on Underside of Buried Structure. The conclusions 
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reached in reference [52] indicate 
that the overpressures exerted on the 
iinderside of buried structures with 
an integral floor are of the same 
order of magnitude as the overpres- 
sures on the roof. Therefore, it is 
recommended ttiat in both the regular 
and Mach reflection regions the bot- 
tom surface of buried structures ^;ltli 
an integral floor be designed to 
support the identical overpressures applied to the roof of the structure as 
given in paragraphs 3- 21a and b. 

3-22 PROCEDURE FOR COMPUTATION OF LOADS OE BUEIED STRCJCTURES. Repeat 
Steps 1, 2, and 3 of paragraph 3-09 . 

Step 4. Determine U , knowing P , from figure 3.9. 

O SO 

step 5* Determine t, . 

d 

Step 6. Determlofi P^ from table 3.I or figure 3.22. 

Step 7. Plot the curves of overpressure vs time as illustrated in 
figures 3.71 to 3.76. 





VS time for buried structure in region of 
Mach reflection 


RADIATION 

3-23 NUCLEAR RADIATIOE PHENOMENA. A characteristic of an atomic detona- 
tion is the emission of nuclear radiation. This represents approximately 
15 per cent of the total energy of a typical air burst. It is expressed 
in terms of initial radiation and residual radiation. 

Initial radiation is defined as the nuclear radiation which is de- 
rived directly from the initial fission eind fusion reactions of the detona- 
tion. It is delivered within approximately the first half minute following 
detonation and its significant effects are confined to a radius of a few 
miles from the point of detonation. 

Residual radiation applies to the radiation emitted after the first 
minute by the fission products, unfissioned residues, and to a limited ex- 
tent by the bomb case fragments and materials such as dust in the air in 
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which radioactivity has heen induced. Residual radioactivity, in. the ag- 
gregate, is of an enduring nature; however, its intensity decreases by a 
process of natural radiation decay. The residual radiation produced by an 
air burst is usually dissipated in the upper atmosphere and does not there- 
fore constitute a hazard on the ground. However, when a nuclear weapon is 
detonated at or near the ground level large amounts of the surface mate- 
rials, such as earth or water, are drawn up into the cloud. Radioactive 
isotopes formed as described above become associated in various ways with 
these materials which, being considerably heavier than those resulting from 
an air burst, will fall back to the earth's surface in the local area; the 
heavier particles falling nearer to the point of detonation and the lighter 
materials settling out progressively farther downwind. The material so de- 
posited is known as "fallout" and the areas on which it falls are said to 
be radiologically contaminated. With megaton yield weapons, falD.out may 
extend several hundred miles downwind and cover five to six thousand square 
miles or more. 

I Fallout patterns are usually depicted as elongated and cigar shaped. 

Such a pattern would only be found under Idealized conditions which seldom, 
if ever, would occur. The existence of cco^Jlicated wind structures, as 
well as variations of these structures in time and space, may cause extreme 
distortions in the shape of the fallout pattern. In addition, intensities 
within the pattern may be extremely Irregular. 

Nuclear radiations consist of neutrons, gamma rays, beta particles, 
and alpha particles. The neutrons, which are subatomic particles of 
neutral charge, are released in the fission and fusion reactions and are 
produced essentially in the first half second after detonation. Gamma 
rays, which are high-energy electromagnetic radiations (like X rays), are 
emitted from the fireball as initial radiation and from fission products 
and from the capture of neutrons in bomb fragments and other materials as 
residual radiation. Beta particles, composed of high-speed electrons 
emitted by the radioactive fission products and alpha particles. Identified 
as the nuclei of helixim atoms, originate in the unfissioned residues of 
plutonium or uranium. Both the neutrons and gamma rays have a long range 
in air arid are very penetrating. Substantial thicknesses of materials of 
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high density are required to reduce their intensity to harmless levels. 

In general, initial gamma radiation is more Important than neutron radia- 
tion for large weapons. However, for weapons of small physical size, 
neutron and gamma radiation may he of equal significance, depending upon 
the characteristics of the shielding provided. Beta particles have a very 
short range in air and relatively little penetration capability. Beta 
emitters can damage human tissue when taken into the body, and in addition 
can cause serious bums if they come in contact with the skin. Alpha 
particles have an even shorter range in air and are stopped by almost any 
type of barrier. They are a hazard only if their emitters are taken into 
the body by inhalation, ingestion or, under limited circumstances, throu^ 
breaks in the skin. 

G amma radiation dosage is measured in terms of a tmit called the 
roentgen r which is a standard measure of the ionization caxised by garmnn 
rays in their passage through matter and hence of the Injury which is 
caused to a body of living organisms. Initial gamma radiation is a func- 
tion of weapon yield, air density, and the distance of the burst from the 
structure. Figure 3.77 gives the intensity of initial gamma radiation at 
any distance frcaa ground zero for a surface burst. 

Neutron radiation dosage, for the purpose of assessing radiation 
hazard, is measured in rem (roentgen equivalent - mammal) which is the 
amount of energy absorbed in mamm alian tissue which is biologically 
equivalent in mammals to one roentgen of gamma or X rays. Neutron radia- 
tion is a function of the weapon yield, weapon design, air density, and the 
distance of the burst from the target. Figure 3.79a gives the iutensity of 
neutron radiation at any distance from ground zero for surface burst of 
fission weapons. Figure 3*79h provides similar data for fusion weapons. 

Nuclear radiation can produce a variety of harmful effects in living 
tissue which may be acute or delayed depending upon the total amount of 
radiation absorbed and the period of time within which it is received. 

When radiation is received within a short time (l or 2 days), the effect is 
considered to be acute and, under most conditions, is essentially inde- 
pendent of dose rate. When dosage is received over a long period of time, 
either continuously or in repeated Increments, partial recovery takes 
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Figure 3.77. Initial gamma radiation, surface burst, 0.9 air density 
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place. Under these conditions of exposure, larger total doses csin be 
tolerated insofar as early effects are concerned. 

Ultimate body damage resiilting from nuclear radiation is a summation 
of the several separate radiation effects involved. Thus, roentgens of 
gamma radiation received must be added to the roentgen equivalent dosage 
for man (rem) of neutron radiation received in order to obtain the total 
dose. The effects of different amoimts of radiation received are shown in 
table 3 •4. 

A particular effect is associated with a range of doses rather than a 
specific dose. This does not necessarily mean that any dose within a given 
range will always result in the effect stated but rather that the dose re- 
quired to cause the effect will fall somewhere within this range. 


Table 3.4. Acute Effects of Whole Body Penetrating on Human Beings 


Dose in 

1 Week, 
rem 

Effect 

150 

No acute effects, possible serious long- tern hazard 

150-250 

Nausea and vomiting within 24 hours, normal incapacitation 
after 2 days 

250-350 

Nausea and vomiting in under 4 hours. Some mortality will 
occur 2 to 4 weeks. Symptom- free period 48 hours to 2 weeks 

350-600 

Nausea and vomiting under 2 hours. Mortality certain in 2 to 

4 weeks. Incapacitation prolonged 

600 

Nausea and vomiting almost immediately. Mortality in 1 week 


In the absence of directives which stipulate otherwise, it is recom- 
mended that protection provided be sufficient to reduce the combined ini- 
tial gamma and neutron dose to not more than 50 rems. This will allow for 
possible later additional dosage received from fallout, within structures 
or outside during decontamination operations, or traversal through con- 
taminated areas. 

3-24 MUCLEAE EADIATIOE SHIELDING. An lnportant consideration in the de- 
sign of a protective structure is the provision of adequate shielding 
against the effects of nuclear radiation. In determining shielding 
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requirements, separate consideration must be given to initial radiation and 
residual radiation. However, ultimate body damage resulting from nuclear 
radiation is a summation of all of the separate radiation effects involved; 
namely, initial gamma, neutron, and residual gamma radiation. 

a. Initial Radiation Shielding. The effectiveness of initial radia- 
tion shielding is a function of the composition and geometry of the shield- 
ing material, the energy distribution of the radiation at the target, the 
distance from the detonation source, the angle of Incidence of the radia- 
tion, and in the case of neutrons, the configuration of the weapon. Ini- 
tial gamma radiation transmission factors for concrete and earth are given 
in figure 3»78* Similar factors for neutrons are given in figure 3*80. 

The curves in these charts are based on the assumption that the radiation 
is essentially point so\xrce and is peipendicular to the slab of shielding 
material. The transmission factor is the ratio of the radiation dose re- 
ceived behind the shielding material to the dose which would have been 
received in the absence of the shield. 

The factors which detennine the effectiveness of an initial gamma 
shield are density and mass. The transmission factors of materials of 
known density other than concrete and earth may be estimated by interpola- 
tion on a density basis, between the curves shown in figure 3* 78* 

Neutron attenuation represents a more complex problem in which 
several different phenomena are involved. First, the fast neutrons must 
be slowed down into the moderately fast range; the moderately fast neutrons 
must then be slowed down into the slow or thermal range; following which 
the thermal neutrons must be absorbed. This absorption process, referred 
to as radiative capture, is accompanied by the emission of relatively hl^ 
energy gamma rays which must be absorbed by the shielding medium. Figure 
3 -80 provides a rough approximation of the neutron attenuation capabilities 
af concrete and earth. More reliable data are not presently available but 
studies in this area are in process. Consequently, neutron attenuation 
data should be checked whenever neutron shielding criteria for si)ecific 
projects are required. 

The procedure for calculating the reqiiired shielding thickness for 
initial radiation is illustrated by the following example: 
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Figure 3.78. Attenuation of initial gamma radiation 
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PjflKrs 3,79b. f-Jeutron radiation dose, surface burst, fusion weapon, 0.9 air density. 
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figure 3.80. Attenuation of neutrons 
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Prollem; A persomel structure is to be designed that will resist a 
50-psi incident oveipressure. A 12- in. -thick reinforced concrete roof will 
meet the stiuctural requirement. Determine the earth cover needed to re- 
duce the initial gw-irmia and neutron radiation to an acceptable level of 
50 rem. Consider a 100-KT surface burst at a distance that would produce 
the design overpressure. 

The weapon effects data are: 

100-KT Surface Burst 

50-psi overpressure - distance from ground zero =0.4 mile 
(figure 3 *121) 

0.4-mile distance - gamma = 70,000 r (figure 3*77) 

0.4-mile distance - neutrons = 100,000 rem (figure 3 •79a) 

Procedure: From figures 3*78 and 3*80, the transmission value in 
feet, of concrete and trial thicknesses for the earth cover are found to 


be as follows; 


Gamma 

Neutrons 


12 in. of 4 ft of 5 ft of 
Concrete Earth Earth 


0.01 

0.003 


0.003 

0.0007 


133 r 


Try 4 ft of earth cover. The amo\mt of radiation transmitted is: 

Gamma: 12 in. of concrete, 70^000 X O.I9 = 13,300 r 

Neutrons: 12 in. of concrete, 100,000 rem x 0.12 = 12,000 rem 

Gamma: 4 ft of earth, 13,300 X 0.01 = 133 r 

Neutrons; 4 ft of earth, 12,000 X 0.003 = ^ rem 

Total transmitted I69 rem 

This is greater than the maximum acceptable level of 50 rem. Try 5 ft of 
earth cover. 

Gamma: 5 ft of earth, 13,300 X O.OO3 = 39 r 

Neutrons: 5 ft of earth, 12,000 rem X O.OOO7 = 8 rem 

Total transmitted 47 rem 

Since this is less than the maximum allowable of 50 rem, 5 ft of earth 
cover is satisfactory. 

b. Fallout Radiation Shielding. The problem of shielding from fall- 
out radiation is different from that of direct radiation. The main differ- 
ences a]re: 
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(1) Eadiatiou frcm fallout particles persists for a long period of 
time "but at a delaying rate. 

(2) The radioactive particles are spread over large areas. 

(3) The radiation received at a point on the grotind comes partly 
from the plane area source and partly from scattered radiations or 
"skyshine. " 

(4) The energy of the radiation from fallout is much lower than the 
Initial gamma radiation. 

Due to the above characteristics of fallout radiation, the calcula- 
tion of shielding thicknesses is more congilex than for direct radiation and 
must take into account the time factor and the hel^t, shape, and size of 
the protected structure together with the location of the point of interest 
relative to the surroimding structure. 

Data and design criteria are contained in references [360 and [49] . A 
detailed method for evaluating the protection afforded by existing buildings 
against fallout radiation is given in reference [50]* This method can be 
adapted for use in the design of shielding for protected areas in new 
structures . 

3--25 TBEEMAL EADIATIOIf. a. Phenomena. The tenqperatures in the fire bedl 
of an atomic bomb detonation are very high (approximately 330 billion calo- 
ries i)er kiloton) , resulting in a large proportion of the energy being 
emitted as thermal radiation. The theimal energy reaching any partlciilar 
point is a function of the distance from the detonation, the yield of the 
weapons, and the scattering and absoiption caused by the atmosphere. The 
aanoxmt of atmospheric attenuation is a function of the visibility, or the 
horizontal distance at which large dark objects can be seen against the sky 
at the horizon. The principal, characteristics of theimal radiation are 
that it; 

(1) Travels with the speed of light ( 186,000 miles per second). 

(2) Travels in a straight line. 

(3) Has very little penetrating power. 

(4) Can be easily absorbed or attenuated. 

(5) Can be scattered. 

(6) Can be reflected. 
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In view of its speed of travel, the time of arrival at a target is altaost 
instantaneous with respect to the detonation. The integrated total thermal 
energy delivered per unit area is plotted in figures 3*81 and 3*82 as a 
function of distance from the point of burst for a clear atmosphere 
(visibility 10 miles) . Intensity of thermal radiation falling on any sur- 
face is reduced in proportion to the cosine of the angle that the surface 
makes with one which would be perpendicular to the incidence of radiation. 
Thus, the intensity on a surface 45 degrees from the peipendicular would 
be about JO per cent of that on a comparable perpendicular surface. 

b. Scaling. If two atomic bombs with total energy yields of 
and Wg are exploded, the proportion of the total energy emitted as 
thermal radiation is approximately the same for each weapon, therefore, 
the thermal radiation energy is proportional to the yield. At a given dis- 
tance from the detonation the total amount of radiant energy on a unit area 
is also proportional to the yield and is shown by the expression 

= Wg/W^ (3.36) 

where and ''^he radiant energy on a unit area for the two 

weapons. 

c. Effect on Materials. A surface exposed to thermal radiation will 
absorb part and reflect part of this radiation. The amount of the absorp- 
tion depends on the material and especially upon its color, as the darker 
colored materials will absorb a much larger proportion than light colored 
materials. The absorption causes the temperature to rise, and damaging 
effects may result. 

The most important physical effects of the high temperatures due to 
the absorption of thermal radiation are, of course, charring or ignition 
of combustible materials and the burning of skin. It is very difficult to 
establish definite conditions under which these effects will or will not 
occur, and their extent; however, for convenience, it is found necessary 
to use total energy per unit area criterion. 

Table 3*5 gives data relative to critical heat energies which produce 
various effects on construction material, surrounding foliage, and exposed 
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Table 3,5>, Critical Thermal Energies, Construction Materials 


Material 

Type 

Damage 

Energy, 
cal/sq cm 




100 KT 

10 MT 

Glass 

Glass reinforced 

Ordinary window 

Bakellte, cellu- 

Highly resistant to heat 
but brittle and will 
break at low blast pres- 
sures (see reinforced 
glass) 

Surface melts and darkens 

TO 

126 

with plastic 
( shatterproof) 

Plastics 

lose acetate. In- 
cite, Plexiglas 

USAF window 

Dense smoking and flames 

430 

750 

Sand 

1/2 in. thick 

Siliceous 

Explodes (popcorns) 

19 

35 

Sandbags 

Cotton canvas 

Chars, bags fail 

18 

32 

Fiberboard 

filled 

Corrugated and/or 

Flames during explosion 

10 

18 

Plywood 

laminated 

Douglas fir 

Flames during explosion 

16 

20 

Wood 

Unpainted planks. 

Deep charring 

18 

32 

Siding steel, 

Douglas fir 

28 gage + 

Not affected 

107 

__ 

corrugated 

Siding- galbe st os 

22 gage 

Flames during explosion 

l 4 

28 

Electric wire. 

I2-2-6OOV 

Flames during explosion 

13 

26 

rubber insulated 

BX armored wire 

600v 

Not affected up to 

107 

16 

Romax nonmetallic 

I9-2-6OOV 

Charring 

8 

Roll roofing 

Mineral surface 

Flames during explosion 

4 o 

71 

Roll roofing 

Smooth surface 

Flames during explosion 

16 

20 

Grass 

Dry 

Flames during explosion 

7 

16 

Leaves 

Dry fallen 

Flames d\iring explosion 

5 

10 
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elements of -structures. It -will "be noted that more heat energy is required 
from the larger weapon than the smaller one to produce the same effect. 

This is due to the length of the delivery time^ which means that in order 
to produce the same thermal effect in a given material, the total amount of 
thermal energy (per unit area) received must be larger for a nuclear explo- 
sion of high yield than for one of lower yield, because the total energy is 
delivered over a longer i>eriod of time. 

Although a serious hazard, injury and damage from theimal radiation 
can be more easily avoided than can some of the other dwtiaging effects of 
atomic weapons. Shelter behind any object which will prevent direct ex- 
posxire to rays from the detonation is sufficient protection. A wall or any 
other screening object, even clothing, will provide a shield from thermal 
radiation. 


CRATERING 


3“ 26 CRATERING PHENOMENA. When an atomic bomb is detonated near the sur- 
face of the ground, a crater is formed in the ground under the explosion 
as shown in figure 3* 83* Two zones of disturbance in the soil aroimd the 
crater are distinguishable. One, the rupture zone wherein the soil has 



Figure 3.83. Typical crater 
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) been violently ruptured, the other, the plastic zone wherein the soil has 
' been permanently deformed but without visible rupture. If the material 
consists of rock there will be a rupture zone but little or no plastic 
zone. Some of the material thrown out of the crater falls back to form 
a lip around the edge. The height of the lip and the extent of the rapture 
and plastic zones are functions of the crater diameter and depths^ as indi- 
cated in figure 3 •83* 

The limits of the rupture zone represent the closest distance to a 
ciatering explosion at which it is feasible to construct a resistant struc- 
I tui^. Due to the detrimental effects of the material thrown out of the 
crater, the higji ground accelerations and displacements, and the high In- 
tensity of the blast pressures and radioactivity, it is not considered 
practicable to design protected structures to resist the effects within 
the radius of the crater lip or within approximately twice the crater 
radius. At distances greater than this, the pressure exerted on buried 
structures by the directly transmitted ground shock is less than the ea 
pressure caused by the air blast wave. Althou^ the directly transmitted 
ground, shock is Initially of higher magnitude than the air blast wave, it 

Is attenuated at a greater rate than the air blast. 

The depth aod dlaa.etar of the crater depend upon the characteristics 

of the soil and upon the energy yield of the explosim as given hy 
figure 3* 8^* 
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Figure 3 . 84 . Crater diameter and depth for a surface hurst. ('Apply factors given above for different soil types) 
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APPENDIX A 

TRANSONIC DRAG PRESSURES 

A-01 INTRODUCTION. This appendix is concerned with the calculation of 
nuclear blast pressures on aboveground earth-covered structures, such as 
shown in figure A.l, which will be applicable to pressure levels above 
25 psig. 

Blast loading ^ ]• 

on structures may i"" ^ 

. FLOW^ 1 

be divided Into two ■ ■ I o 

stages, i.e. the | h ^ 

diffraction phase ^ structure 

62' 62' 

and the enveloped r 

phase . When the 

faj Case 1 

pressure at the 

rear of the surface i 

I 16 1-16 -I 

attains the value r — T \ 



1 

_ _ _j 


A 

62' 

(^structure D 
62' 


(a) Case 1 

i 


of the pressure in FLOW ^ 
the blast wave, the 

diffraction process A (^structure ^ 

may be considered 56^^ ^ »■ 

to have terminated, Case 2 

pir\^ subsequently^ 

Figure A.L Aboveground earth-covered structures 

steady state condi- 
tions may be assumed to exist until the pressures have returned to the am- 
blent value prevailing prior to the arrival of the blast wave. The dif- 
fraction phase pressures on the surface of the earth cover may be computed 
by methods set forth in EM 1110-345-413 and EM 1110-345-420; however, little 
information is available concerning the enveloped phase pressures. The en- 
veloped phase pressures consist of the incident pressures plus the dynamic 
pressures. This appendix presents a solution for obtaining the enveloped 
phase pressures on wedge-shaped earth covers similar to those shown in 
figure A.l. The pressures transmitted through the earth to the structure 
may be computed by methods described in EM 1110-345-413 and EM 1110-345-420 


(^structure ^ 

56' 56' 

(b) Case 2 

Figure A.l. Aboveground earth-covered structures 
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or iDy use of a Mohr's Circle Solution. The results are ccrapared with the 
equivalent circle solution given in these ma n uals . 

The flows considered are those "behind the shock waves whose charac- 
teristics are given in table A.l. 

Table A.l, Shock Wave Characteristics, 25 psig to 200 psig 


Shock pressure, psig 

25 

50 

100 

200 

Dynamic pressure, psig 

12.2 

kl 

123 

330 

Shock velocity, ft/sec 

1756 

2216 

2927 

3985 

Sound velocity behind shock, ft/sec 

1309 

i¥j -7 

1681 

2064 

Particle velocity behind shock, ft/sec 

868 

1375 

2082 

3059 

Temperature at shock, deg Kelvin 

379 

467 

636 

968 

Density at shock, p x 10*^ Ib/in.^ 

871 

1161 

1526 

1895 

Mach number behind shock 

0.661^ 

0.951 

1.239 

1.481 

Stagnation pressure behind shock, psla 

53.5 

116 

294 

766 


The flows and the structures are considered two dimensional. This is 
conservative with respect to the pressures on the structure. The effect of 
the structures being finite in the direction perpendicular to the plane of 
figure A.l is to reduce somewhat the pressures near the ends of the struc- 
ture as compared to the pressures computed in the following paragraphs. 

A-02 GEHERAL FEATUEES OF THE FLOW. The general features of the flow de- 
pend on the Mach number of the given conditions behind the original shock 
wave. The Mach niimber is the ratio of the particle velocity to the veloc- 
ity of sound. When all parts of a flow have a Mach number less than one, 
the flow is called subsonic. If the Mach number is everywhere greater than 
one, the flow is called supersonic. If the Mach number is less than one in 
some places and greater than one in other parts of the same flow, it is 
called transonic. All of the combinations of initial flows and the 
geometries given fall in the transonic range for which both experiments 
and theory are, unfortunately, relatively scarce and complex. 

In the subsonic case the usml concept of pressure coefficients Cp 
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to be imiltiplied by the dynamic pressure l/2pV^ are applicable. Coeffi- 
cients given for incompressible flow may be used with less than 1 percent 
error for Mach numbers from zero to 0.1. This covers the range intended 

for most wind load coefficients. From Mach numbers 0.1 to 0.5 the incom- 

1 


pressible pressure coefficients may be used if first atultiplied by ^ . 

This is the Prandtl-Glauert rule. (See reference [a. 2], p 139). For^ ' ^ 
larger Mach nxmibers, more elaborate methods must be used. 

In the supersonic case the flow pattern is illustrated in figure A. 2. 



On hitting the comer at A, the flow forms a shock wave which turns the on- 
coming flow abruptly, so it is parallel to AB. This causes a pressure 
rise. At B, the flow goes through a Prandtl-Meyer expansion fan which 
again turns the flow, this time with a pressure reduction. At C there is 
another expansion fan and then a compression shock wave at D returning the 
flow to its original direction. The changes of pressure, density, tenipera- 
ture , etc . , at each shock and expansion can be conveniently figured using 
figures A. 3 through A. 5 . The pressure is relatively constsint on each of 
the faces AB, BC, and CD, being of course greatest on AB. The shock wave 
at A is called an attached shock wave because it starts right at A. As the 
Mach number decreases toward one, there will be a certain Mach number 
(figure A. 6) at which the attached shock will no longer be possible for the 
given angle 9. Below this Mach number the shock wave is detached as shown 
in figure A . 7 and the flow will be in the transonic range. 

In the transonic flow shown in figure A. 7, the flow is supersonic 
throughout except for the embedded subsonic region EFBAE. The initial 
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Figure A.3. Ratio of supersonic stagnation pressure to free stream pressure vs Mach number 



Figure A.4. Deflection angle required to expand to a given Mach number 
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flow is assumed stqpersonlc, 
l.e. M > 1. (Reference 

[A.T].) 

As the initial Mach 
number decreases to 1, the 
detached shock in figure 
A. 7 forms further ahead of 
point A. When the Mach 
number falls below 1, the 
typical situation in the 
transonic range is as 
shown in figure A. 8 with 
one or more supersonic xeglohs embedded in the subsonic flow. The condi- 
tions shown apply to Mach numbers not too close to one. As the Mach num- 
ber increases and approaches one, the supersonic regions grow and coa- 
lesce, finally covering the structure from B to D. In this condition 
there would be a shock wave starting from D as in figure A.J. 

The flows shown in figures A. 7 and A. 8 have certain features of 
interest in common. In both cases the pressure at A is equal to the stag- 
nation pressure of the subsonic flow upstream of A. At B the local Mach 
number is always unity at the start of the expansion fan. These facts have 


M-Mach Number 

Figure A.6. Maximum deflection angle for which the 
shock wave will remain attached (reference^^) 


F 
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Figure A. 8. Transonic flow (M^ < 1) 


befen establisbed experimentally and theoretically (reference [A.6]). 

Further, it has been shown that the pressure distribution, on the face 
AB is the same for a wide range of transonic flows for a given angle 9 
when the pressure is expressed in terms of the stagnation pressure rather 
than a pressure coefficient (reference [A.53^ P 255) •*■ These facts are the 
basis for the method used below in computing the pressure distribution j^or 
the cases at hand. 

a -03 method for computation of pressures ii^ the transonic range . 

a. Cases In Which The Initial Flow Is Supersonic . The two highest 
pressure shock waves of 100 psi and 200 psi shock pressure yield an initial 
flow which is supersonic as shown in Table A.l. In these cases the condi- 
tions shown in figure A. 8 apply. 

The first computation is to find the stagnation pressure in the sub- 
sonic region. This may be conveniently done using figures A. 3 and A.-5. 

This stagnation pressure at point A, (figure A. 9) is obtained by first 
determining the supersonic downstream pressure, Psi, corresponding to the 
freestream Mach number, M^ = M^, from figure A. 3 (using Y = 1.4), then 
using figure A. 5 to determine the 
pressure at A, Ps2. For points be- 
tween A and B, for the particular 
angle involved (0 = 26.6®) there are 
fortunately both theory and experi- 
ments available giving the ratio of Figure A. 9. Definition of x and c 
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the pressure at any point x to the stagnation pressure which occurs at A. 
These ratios ane plotted in figure 16 of reference A. 6 and are summarized 
in Table A.2 


Table A»2, 


Pressure ^ ^ 26.6 ° 

Pressure at A, ps 





It will be noted that the experiments and theory in most of the refer' 
ehces deal with wedges as shown in figure A.IO. The structures considered 

herein are effectively one-half of the 



Figure A.IO. Typical wedge profile 
not available for the faces BC and 


wedges tested. 

The test results and theory avail- 
able are concerned primarily with the 
face AB which is the most heavily loaded 
region and the relatively precise data 
given in Table A.2 is applicable. Un- 
fortunately, equally precise data is 
). On the basis of examination of test 


results and of the underlying theory the following procedure is suggested as 


conservative auid adequate, though lacking the same precision as the loads on 


AB. 

The pressure at B, on the BC side of the corner may be computed by 
assuming there is a Prandtl-Meyer expansion fan starting from Mach number 
one on the downstream side of the corner as the net effect of an over- 
expansion and compensating compression at the corner. This pressure may be 
obtained by determining the Mach number, on the BC side of the corner 
corresponding to the deflection angle Q = 26.6° from figure A. 4 and then 
entering figure A. 3 with this value to obtain the pressure, P^. Along the 
face BC there may be some weak or even quite strong shock waves for which 
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exact data is not available, but their effect will be to increase, the 
Ijressure. A reasonable estimate for the cases being considered is tha.t 
the flow is returned to the pressure of the original supersonic flow up- 
stream of E by the time the corner C is reached. In the absence of exact 
data a linear distribution of pressure between B and C is assumed. 

At C, the pressure on the downstreaun side may be estimated by 
assuming a Prandtl-Meyer expansion fan starting from the free stream 
Mach number. Along the face BD, it is suggested that the pressure be 
assumed uniform because this region is shielded from the influences that 
can cause pressure rise along BC. The pressure on CD may be obtained 
from figure A. 4 and A. 5a as follows: (1) enter figure A. 4 with free 
stream Mach number to obtain the Prandtl-Meyer angle, 0, (angle tha^u 
which a supersonic stream is turned to expand from M = 1 to M?’!), (2) 
add 0 to the deflection angle ( 26 , 6 ® for our case) to obtadn the ratio 
p/Pt where p is the static pressure on face CD and Pt is the total pressure 
downstream of C. Since there are only isintropic processes between A and 
D, Pt is the same as at A, (3) compute p from the ratio p/Pt. 

b. Cases In Which The Initial Flow is Subsonic . The shock waves of 
25 and 50 psi in Table A.l have original Mach numbers less than one. In 
these cases, the pressure at A is the stagnation pressure given in Table 
A.l and is obtained from figure A. 3 only by entering the figure with the 
free stream Mach number to obtain the pressure ratio. 

The coefficients given in Table A. 2 apply for the region AB, using 
the stagnation pressure as given in Table A. 2 for the angle 0 = 26.6® 

The pressure on the downstream side of B is computed assuming an 
expansion fan from Mach number one and may be obtained from figure A. 3 
and A.4 for 0 = 26.6® as for the supersonic case. A linear pressure rise 
is assumed from B to C, reaching the free stream pressure at C. 

The downstream pressure at C is estimated according to their being 
a separation of the subsonic flow at the corner and with th« help of 
experimental data (reference A.9«) Simple inviscid theory would indi- 
cate that the flow expands to zero pressure at the corner and thereafter 
resumes free stream velocity, Hov/ever, since for these Mach numbers of 
approach to an edge of such a large angle as that at C, the flow separates, 
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viscous effects can not be ignored (reference A. 9). The pertinent data 
in reference A. 9 (page 28 for station 12, l6) applies to a 30 ° flap de- 
flection at an angle if attack of the flap of -9“. It would be desira- 
ble to have data for a leading edge flap deflection of 26.6® and an 
angle of attack of the flap of 0® but the above data is the closest 
available at this time. The pressure coefficient, P = (p-PQ)/q indi- 
cated by the above test data is between -0.6 and -0,7, where p = the 
static pressure (psia), P = free stream pressure (psia) and q = the 
dynamic pressure. From another standpoint, if we can apply data from 
a wedge to represent flow on a half wedge following a straight section, 
as we did on the forward face AB, then we can apply data from a sharp 
edged plate at an angle of attack of 26.6® to represent the flow, around 
a sharp corner of 26.6®,- with about the same validity. The data on 
pages l8 and 19 (zero flap deflection) of reference A,9 show that the 
flow is separated and P is constant over the entire chord for I50, so 
it surely is also for «= 26.6®. The value of P is about 0.6, and this 
is used for both the 25 and 50 psi cases in this report, 

A-04 NUMEIBICAL BESULTS FOB THE GIVEN CASES. The absolute pressures in 
pounds per square inch (absolute) are shown in figure A. 11 for the given 
shock waves of 25, 50 , 100, and 200 psig, respectively. The exact varia- 
tion of the pressure along AB in each case may be found by applying the 
coefficients of table A. 2 to the pressure given for point A (figure A. 9 ) 
in 6ach case. For comparison, the pressure distribution obtained using the 
equivalent circle method described in EM 1110-3^5“^13 EM 1110-3^5-^20 
is shown for the 25-psig level, and the 200-pslg level in figures A. 12 and 
A. 13. The values given in figure A. 11 are based upon a "flat top" shock 
but the preceding method is also easily applied to shocks which vary with 
time. 

The pressures given are intended to apply to both of the structures 
shown in figure A.l(a) and figure A.l(b). The difference in the length of 
BC in the two cases will, in fact, make little difference to the pressure 
distribution and the available information does not provide any basis for 
computing these variations. It would be desirable to have a solution for 
flow over a double-wedge profile of 53° nose angle, but we have not been 
able to locate any data for this range. 
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PRESSURES, PSI ABSOLUTE, FOR 25PSIG SHOCK WAVE 



PRESSURES, PSI ABSOLUTE, FOR 50PSIG SHOCK WAVE 



PRESSURES, PSI ABSOLUTE, FOR lOOPSIG SHOCK WAVE 



PRESSURES, PSI ABSOLUTE, FOR 200PSIG SHOCK WAVE 


Figure A. 11. Steady state pressures for the structures in figure A.l 
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WEDGE METHOD 


Flo# 39.7 


■ structuhc. 



\\\ i 


PRESS. IN FRONT 
OF PT © 8 TO REAR 
OF PT@ NOT REQUIRED 
TO OBTAIN PRESS. ON 
STRUCTURE 


EQUIVALENT CIRCLE METHOD (EM 1110- 345- 413 AND EM 1110-345-420) 
note: all pressures are psi absolute 


Figure A.12. Comparison of wedge method with equivalent circle method, pressure - 25 psig 
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Flow 214.7 




EQUIVALENT CIRCLE METHOD (EM IIIO-345-4I3 AND EM 1110-345-420) 
NOTE all pressures ARE PSI ABSOLUTE 


Figure A. 13. Comparison of wedge method with equivalent circle method, pressure = 200 psig 
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■A-05 KEMASKS COHCERITIIirG DOME OR COIflCAL- SHAPED STHJCTUP'P^P - The above 
analysis has been concerned with structures which are long enough con^iared 
to the dimension in the direction of flow to be considered two dimensional. 
A dome or conical- shaped structure does not fall in this category and the 
coefficients suggested above should not be used for such cases, although 
this would certainly be conservative as far as the heavily loaded windward 
side is concerned. 

Field measurements on full scale structures of this type were made 
during the 1957 nuclear test series conducted at the Nevada Test Site, 
but were not available at the time of this report. 
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